Netherlands 
organization  for 
applied  scientific 
research 

TNO-report 


’•yoit  no-  copy  no.  _ 

FEL-91-B112  y 


AD-A245  601 

illliiill 


TNO  Physics  and  Electronics 
Laboratory 

RO.  80x96864 
2509  JG  The  Hague 
OvKle  Waalsdorperweg  63 
The  Hague.  The  Netherlands 

Fax  +31  70  32609  61 

^  Phone +3 1  7032642  21 

.Radio  chanoel  measumiieDts  usiog  a  tarosAasd 
PMudo-Noue  tignal  (ffleasareaieiit  aet-op  and 
procesaing  of  the  lesuhs) 


Nothing  from  this  issue  may  be  reproduced 
and/or  puCHished  by  print,  photoprint. 
microWm  or  any  other  means  without 
previous  written  consent  from  TNO, 
Submitting  the  report  for  inspection  to 
parties  directty  interested  is  permitted. 


auVior(a): 

Jr.  J.AM.  Vriens 


In  case  this  report  was  dratted  under 
instruction,  the  rights  and  obligations 
of  contracting  parties  are  subject  to  either 
the  'Startdard  Conditions  for  Research 
Instructions  given  to  TNO"  or  the  relevant 
agreement  concluded  between  the  contracting 
parlies  on  account  of  the  research  object 
involved 

irTNO 


Ir.  GJ.M.  Janssen 


June  1991 


danfficasan 

«• 

abdnct 

npoct  «ixt 


MC 


:  unclassified 
unclassified 
:  unclassified 
;  unclassified 


92-02827 

■iiiimi 


no.a(cautM  ..  25 

ocL  RI» dMflnkn  Uip 

•fipmtaum  ;  3 


M  Monvudon  wtikl\  la  daaillad  UBeanlng  la 

Oiaoh  rauiiilpna  «WI  ba  auilau  by  in 

***  — "  "W  aa  daaaHad  miotinWun  el 
ow"Won«i9  vik»  kt  Ma  Don  eeurky-  No  yat  el 
C*  Nknmaen  «■  ba  «MlDaad  ID  aiy  paly. 


CD 

CO 

-f 


oz 

o 

CM 

CM  j 

A  i 


TNOrapoit 


Page 

2 


raportno. 

tWe 


aulhor(s) 

Intittuto 


NORO  no. 
no.  In  pow  ai 

RoMareh  suparviaed  by: 
naiaaroh  earrlad  out  by : 


FEL-91-B112 

Radio  channel  measurements  using  a  broadband  Piendo-Noise  signal  (messatemem  Kt-up 
and  processing  of  the  results) 


b.  Vriens,  b.  GJ.M.  Janssen 
TNO  Physics  and  Electronics  Laboratoiy 

June  1991 

711.2 


Accesion  For  ^ 

NTIS  CRA&I 
OTIC  TAB 

U  a.ii'.Oiji  CL-d 
Jl;Stil;CatlO:l 

1/ 

By . 

D;.  t  pbjtic  ./ 

Diit 

A') 

/* 

{ 

b)  a  multipath  communicatioiis  envbonment,  which  is  often  found  in  uiban  areas,  bansmission  of 
infonnation  is  difBcub  due  to  self-inteifeience.  This  self-interference  may  lesuU  in  signal  sqtpression  and 
intersymbol  interference.  Knowledge  of  die  parameters  that  characterize  the  radio  channel,  can  be  used  to 
feed  an  adaptive  equalizer  or  frequency  management  system,  in  order  to  minimiM  diese  disturbing  effects. 


bi  this  report  a  measurement  system,  to  measure  the  characteristics  of  radio  paths  in  a  multqiath  channel,  is 
described  and  measurement  results  ate  ghren. 

The  measurement  system  uses  Pseudo-Noise  Spread  Spectrum  modulatko  and  operates  in  die  VHF-band 
(30  -  70  MHz).  A  cania,  BPSK  modulated  widi  die  Pseudo-Noise  signal  at  a  rate  of  1  Mch^/s,  is  used  to 
sound  the  radio  channeL 

The  characteristics  of  the  different  radio  padis;  relative  amplitude,  RF-phase  and  delay-time,  are 
from  die  coberendy  received  signaL 

The  accuracy  of  die  measurement  system  has  been  determined  and  is:  2  dB  for  die  umpiinnU.,  50  for  the  KP- 
phase  and  20  ns  for  die  debqr-time. 
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With  the  system  measmements  have  been  perfonned  in  a  contioUed  multqtath  envinmment  in  die  labwatoiy 
and  in  the  field  in  different  types  of  teirain.  The  contioUed  multqiadi  results  are  presented  in  diis  tq;Kitt,  die 
analysis  of  the  measurements  performed  in  different  types  of  tetrain  will  be  presented  in  [Vriens,  1991]. 
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SAMENVATTING  (ONGERUBRICEERD) 

In  een  omgeving  met  meerweg  radiqnopagatie,  zoals  vaak  in  siedelijke  omgevingen  vooikomt,  wordt  bet 
veizenden  van  infonnatie  bemoeilijkt  door  zelf-interferentie.  Deze  zelf-interfoentie  Van  lesulterea  in 
signaalondefdrukidng  en  intersymboolinteiferentie.  Indien  de  paiameten  die  bet  nditdtanaa)  karakterisetcn 
bekend  zijn,  kan  deze  infonnatie  wcnden  grixoikt  voor  bet  sturen  van  een  adsptief  egalisadefilter  of  een 
frequentiemanageinent  systeem  om  zo  de  negatieve  eSecten  te  mmitiMiiiinmn 

In  dit  rapport  wordt  een  meetsysteem  beschteven  waannee  de  kaiakteristidcen  van  ladit^Mden  in  een 
meerwegomgeving  kunneu  worden  gemeten.  Tevens  worden  mwrf  gegeven. 

Het  meetsysteem  maakt  gebiuik  van  "Pseudo-Noise  Spread  Spectrum”  modnlatie  en  weikt  in  de  VHF-bHid 
(30  -  70  MHz).  Een  draaggolf.  BPSK  gemodnleenl  met  bet  ”Piendo-Noise"-signaal  (met  een  chqt-rate  van 
WQidt  als  zendsignaal  gefaraikt 

De  karakteristiekea  van  de  verschnieade  radiiqraden:  relatieve  anylitnde,  RF-tee  en  venragingstijd,  fcimiinn 
ttit  bet  coberent  ontvangen  signaal  worden  berdtand. 


De  nauadceori^ieid  van  het  meetsysteem  U  bqpaald,  en  is:  2  dB  voor  de  ampHttvy,  50  voor  de  RF-fine  en 
20  ns  voor  de  vomgingstijd. 


TNOraport 


Pl^ 

5 


Met  het  systeem  zijn  metiiigen  uitgevoerd  «an  een  geconaoleenle  meerweg  omgeving  in  het  labantofinm  ea 
in  het  veld  in  verschillende  typen  tenein.  De  re*nlt#ten  van  de  geconlToleenie  meerweg  omgeving  worden  in 
dit  rtf>poit  gegeven.  De  analyse  van  de  metingea  oitgevoerd  in  verschillende  typen  tenrein,  zuUen  in  [Vriens, 
1991]  worden  gepresenteerd. 
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1  INTRODUCTION 

In  todays  conununicatioiis  the  rote  of  mobile  commiinicatiQns  is  increasing  nqndfy.  In  a  mobile 
/vwTiitiiiiiifiirinn«  environment  however,  the  transmission  channels  are  not  optimum;  operation  has 
to  be  petfmned  in  a  more  or  less  severe  multipath  environment  and  among  odier  interfering 
signals. 

lYopagation  of  radio  signals  between  a  tmtsnuOa  and  receiver  strongly  depeitds  on  die  gemnetry 
of  the  radio  channel  enviraunent  When  nmltipath  propagation  occurs  (which  is  caused  by 
different  mechanisms  as:  reflection,  diftactiiw  and  scattering),  the  desired  signal  in  general 
reaches  the  receiver  via  more  tiian  one  propagation  padi.  Due  to  the  occurring  differences  in 
mnplitiiHe,  delay-time  and  RF-phase  of  the  different  paths,  the  received  signal  will  be  strongly 
distorted  or  extinguish  totally  at  certain  frequencies,  t^ch  makes  communication  difficult  or 
iiiqtossible. 

When  the  multipath  characteristics  are  known,  this  information  can  be  used  to  counteract  these 
negative  effects. 

For  hi^  data  rate  digital  transmission  systems  (e.g.  trunk  connections)  tire  multipath  infrjtmaliiHi 
can  be  used  to  feed  an  adaptive  equalizer. 

For  narrow  band  operation  the  multqratii  information  can  be  used  to  calculate  the  frequency 
response  of  a  imilripath  channel  in  order  to  choose  fee  tqrtimum  frequency.  For  proper  operation 
of  the  adiqitive  equalizer  and  in  order  to  choose  the  optimnm  frequracy,  accurate  channel 
infonnatioa  is  necessary  on  relative  arrqrlitude,  RFffease  and  delay-time  of  tire  occoning 
multqMife  signals. 

In  this  rqrort  a  channel  measurement  tedmiqne  is  presented  wUdi  makes  it  possiMe  to  measure 
tile  conqtiex  amplitudes  and  delay-times  of  tire  radio  arriving  via  different  propagation 

paths.  The  measurement  system  is  designed  to  operate  in  tiie  VHP  (Very  Hig^  Frequency)  band 
from  30  -  70  MHz.  The  radio  channel  is  sounded  by  means  of  a  btoaifi>and  Fsendo-Noiae  (FN) 
signal  At  tiie  receiver  side,  die  conqriex  ampHtode  and  delsy-time  of  die  differeM  paths  can  be 
determined  by  means  of  soAware  imptemesited  correlation  f**^^"** 
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The  measurements  to  be  performed  with  die  presented  technique  serve  die  following  goals: 

To  detetmine  the  statistical  characteristics  of  the  multipath  radio  channels  (delay-time, 
amplitude)  for  different  terrain  types; 

To  investigate  die  possibility  of  using  multqiadi  diaracteristics  to  feed  an  adtpdve  equalim 
in  high  rate  data  links; 

To  investigate  the  possibility  of  using  nmltipath  characteristics  for  evaluating  die  narrow¬ 
band  communicadon  possibilities  over  a  wide  frequency  band. 

In  this  report  die  following  items  will  be  discussed: 

In  chapter  2  die  model  of  die  multqiath  radio  channel,  as  used  in  diis  report,  is  defined. 

In  chapter  3  different  princ^Ies  of  channel  measurement  techniques  are  discussed  and  the 
broadband  channel  measurement  technique,  as  used  here,  is  described. 

•  In  cluqiter  4  the  Pseudo-Noise  measurement  technique  is  described  in  more  detail. 

In  dupter  3  the  processing  of  the  received  signal,  which  has  been  inplemented  in  software, 
is  given. 

In  chapter  6  the  accuracy  of  this  technique  is  investigated. 

-  In  chtqiter  7  the  analysis  of  the  narrow  band  fiequency  response  of  die  communicatiim  links 
under  consideraticm  is  described. 

In  chapter  8  the  results  of  tests  in  a  siinulated  multipath  environment,  ate  discussed, 
iqnally  conclusions  and  recommendations  are  given. 

Results  of  the  measurements  which  have  been  performed  in  different  types  of  tenain,  will  be 
presented  in  [Vriens,  1991]. 


s 

(! 
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2  DEFINmON  OF  THE  RADIO  CHANNEL 

Very  often,  eqtecially  in  an  tnban  environment,  die  propagation  of  radio  signals  between  2  points 
retrxi  piw  via  different  paths;  a  direct  path  and  one  or  more  indirect  padis.  This  is  called 
imiiripaih  propagation.  Ih  ftgute  2.1  a  multipath  sitnation  is  sketched. 


Hg.  2.1:  A  schematic  view  of  a  moltipalfa  channel 

Different  propagatiaa  effects  diat  may  ph^  a  role  ate:  reflection,  dif&actian,  deby-time,  signal 
attenuation,  etc. 

When  die  radio  diannel  consists  of  only  one  padi,  die  radio  signal  is  only  attenuated  and  delayed. 
When  the  transfer  ftinction  is  fieqnency  independent,  no  distortion  will  occur.  Ifowever,  in  a 
mnltipadi  situatioa  fieqnency  dependent  constructive  and  destructive  intetference  will  occur  due 
to  summation  <rf  different  attennwd  and  delayed  signal  conqionents.  This  may  cause  severe 
signal  distottian  or  even  total  extinction  of  the  signal 
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When  the  multipath  channel  is  known  accurately  over  a  broad  frequency  range,  measures  can  be 
taken  to  inqnove  communication  possibilities.  In  broadband  communication  systems  the 
distorting  multipath  effects  can  be  ctmqtensated  by  using  equalization.  For  nanow  band  ^sterns 
optimum  frequency  selection,  witii  maTimmn  power  transfer  (cmistmctive  interference)  is 
possible. 

In  multipath  propagation  a  number  of  diffierent  mechanism  play  role.  When  all  these  mechanisms 
are  taken  into  account,  the  desct^on  of  die  radio  channel  will  become  very  conq>licated. 
Therefore  a  number  of  assumptiims  will  be  made  for  the  model  of  die  multipath  radio  channel  that 
will  be  used  tfarou^out  this  tepott 

1.  The  multipath  channel  is  a  time  independent  radio  channeL  This  is  a  valid  assunqition  for 
communication  systems  with  ntm-moving  transmitter  and  receiver. 

2.  For  a  single  path,  pnqiagatim  is  assumed  to  be  frequency  indqiendent  This  means  that 
attenuation  and  delay-time  are  assumed  constant  as  function  of  frequency  (no  dispersion). 
This  assunqition  results  in  a  mayor  sinqilification  of  the  model  In  practice,  however, 
attenuation  varies  as  201og(f)  with  frequency.  But  for  bandwiddis  within  10%  of  the  cairier 
frequency  it  is  a  reasonable  assunqition. 

3.  The  change  of  die  RF-fdisse  on  different  propagation  padis  as  functitw  of  frequency,  is  only 
caused  by  die  differences  in  padi  length  (or  delay-time). 

4.  To  limit  the  number  of  padis  in  a  practical  multqiadi  situation,  only  prqpagatitm  paths  widi 
less  than  30  dB  attenuation,  when  conqiared  to  die  strongest  signal,  ate  assumed  to  have  a 
significant  effect  on  die  total  channel  reqwnse. 

Based  on  these  assumptions,  die  multqiadi  diannel  can  be  modelled  as  the  impulse  response  widi 
complex  transfer  function  h(t),  which  is  die  sum  of  the  sub-transfer  functions  belonging  to  every 
single  padi: 


N  IT 

h(t)  -  Sft-Hi)  (2.1) 

Where  Aj  it  the  ampKtnde,  4)'  is  die  RF-pbaae  diift  doe  to  reflections  (in  order  to  distingntdi  it 
from  the  phase  shift  due  to  the  path  delay-time  xj),  and  xj  ibe  delay-time  of  the  padi.  8(t-Xj)  is 
the  Dirac  impulse  fltnction  at  jtistant  Xj.  N  is  the  number  of  letevant  mnttipmh  chanods. 
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When  the  complex  trananitmd  and  received  signals  are  indicated  as  s(t)  and  ift)  reqrectively 
(with  real  parts  8j(t)  tesp.  r^Ct)  and  imagiiwry  parts  SjCt)  resp.  rjCt)),  the  signal  arriving  at  the 
receiver  for  this  type  of  channel  is  given  by  i^t): 


r^tt)  -  Raal{a(t)  *  h(t)) 

K 

«  Real(  Z  Aj,  expCj^j^')  B(t-ti)) 


(2.2) 


which  is  the  sum  of  attenuated,  phase  shifted  and  delayed  versions  of  the  ariginal  transmitted 
signal.  Every  single  path  propagates  the  signal  undistorted. 

If  (2.2)  is  analysed  it  follows  that  die  received  signal  is  completely  determined  by  the  relative 
amplitude-,  RF-phase-  and  delay-time  values  of  the  different  paths  in  the  mult^th  channel. 
Knowledge  of  these  parameters  allows  a  perfect  prediction  of  the  radio  charmeL 
Theoretically  this  is  true,  however,  in  practice  these  parameters  can  only  be  determined  with 
limited  accuracy. 

It  will  be  shown  that  (widi  respect  to  the  delay-time  accuracy)  for  a  good  frequency  response 
prediction  over  a  reasonable  bandwidth,  it  is  sufBcient  to  know  the  relative  delay-time  values  with 
limited  accuracy,  when  the  relative  RF-phase  values  of  the  multqMth  comptments  ate  known 
accurately. 
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3  CHANNEL  MEASUREMENT  TECHNIQUES 

Several  measurement  techniques  exist  to  detennine  the  multipath  characteristics  of  a  radio 
channel.  These  multqiath  characteristics  are  the  number  of  relevant  paths  from  transmitter  to 
receiver  and  their  reqiective  amplitude,  RF-idiase  and  delay-time  differences,  as  discussed  in 
chapter  2.  All  paths  together  form  the  radio  channel  with  impulse  response  h(t): 


N  N 

h(t)  -  X  hi<t)  -  I  AiOxp(j«|>i')S(t-ti)  (2.1) 

i-1  i-1 

In  order  to  measure  the  multipath  characteristics  of  the  radio  channel,  a  number  of  techniques  are 
available.  In  all  these  techniques  the  radio  channel  is  soimded  by  a  carrier  signal,  modulated  with 
a  specially  chosen  waveform,  which  makes  it  possible  to  distinguish  the  contributions  of  the 
different  paths  in  the  total  received  signal.  In  the  following  tiiree  broadband  noodulating 
waveforms  will  be  discussed  briefly: 

Pulse  modulation. 

Frequency  chirp  wave  modulation, 

Pseudo-Noise  modulation, 

3. 1  Pulse  modulation 

With  this  technique  a  pulse  of  very  short  duration,  modulated  on  a  carrier  signal,  is  transmitted. 
The  complex  transmitted  signal  can  be  indicated  as: 

8{t)  -  F(t)exp(j(Dct)  (3.1) 

with: 

<D^  is  the  carrier  frequency 

P(t)  is  the  anqriitude  modulating  pulse  signal;  P(t)  >  U(t)  -  U(t-‘^) 

U(t)  is  the  unit  step  frinction 
is  the  duration  of  the  pulse. 


The  actual  transmitted  signal  is  now: 
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Sj.(t)  »  P(t)COS(<ll>Qt) 

With  (2.2)  it  follows  for  the  received  signal; 


(3.2) 


M 

r^(t)  -  exptj^j^')  Plt-T^)  oxp(  jlDj.(t-ti) )  )  (3.3) 

In  the  ideal  case,  this  pulse  has  an  infiiiitely  short  width  and  Aerefore  occnpies  infinite 
bandwidth.  If  we  assume  that  all  paths  have  infinite  bandwidth,  for  each  padi  a  similar  resptmse  is 
received.  The  received  signal  may  look  like  figure  3.1.  From  this  figure  the  parameter  values  for 
the  relative  an^litude  and  delay-time  can  directly  be  derived. 


In  practice,  however,  aeveial  ptobknis  ariae  to  in^Ienieiit  fiiis  techniqae. 

Rrst,  we  cannot  generate  a  pulse  widtinfinitdy  short  duratkn.  Consequently,  two  padis  may  hsnre 
a  delay-difference  less  than  die  pulse  duration  and  received  pulses  will  mterfere.  fia  figure  3.2  diis 
situation  is  dqiicted  in  case  two  path  reqionses  interfese.  It  takes  more  effort  to  separate  die  padi 
delay-diiies  of  diese  pulses. 
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Fig.  3.2:  Pulse  response  overlr^  due  to  multqMth  inteifeience 

In  pcacticc,  the  pulse  will  have  a  duration  in  die  order  of  1  microsecond.  The  pulse  is  rectangular 
and  modulates  a  earner  at  a  cotain  frequency  in  die  VHF-band.  Dqiending  on  the  exact  pulse 
shape  and  duratimi,  the  pulse  occiqiies  a  bandwidth  of  several  megahertz,  centred  around  the 
carrier  frequency. 

The  second  problem  concerns  the  transmitter.  Since  die  transmitted  pulse  has  a  very  short 
duration,  a  large  amount  of  energy  must  be  concentrated  in  a  very  diort  time,  in  order  to  achieve  a 
good  signal-to-noise  ratio  at  die  receiver.  Therefore,  a  very  large  peak  power  must  be  available  at 
the  transmitter.  This  puts  hi^  demands  on  the  power  amplifier  of  die  transnoitter. 

The  third  problem  widi  this  tedbuique  concerns  the  receiver.  Since  most  of  the  time  onty  noise 
and  no  signal  is  present,  die  receiver  must  have  a  large  dynamic  range  to  process  die  unexpected 
pulse.  IHodiennore,  (he  duration  of  die  pulse  makes  it  inqiossible  to  receive  die  signal  crdierentiy. 
Hence  no  phase  infonnstion  is  detected  and  only  die  energy  of  each  padi  instead  of  the  conqdex 
amplitude  can  be  determined. 


! 
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3.2  Requency  chirp  wave  moduladon 

In  this  technique  flte  transmitted  signal  is  a  caniei  signal,  v^iich  is  poiodically  swq>t  in 
frequency  from  to^jn  to  'The  complex  transmitted  signal  can  be  written  as: 

a(t)  -  +  At/2)t)  (3.4) 

with; 

A  is  the  sweep  rate  (Wnum  -  [rad/s^] 

T  is  the  sweep  duration. 

The  signal  arriving  at  the  receiver  is  now; 


H 

r^lt)  -  Real{^X,^AiOxp(J6i')«xp(j(«Binin  +  A(t-Xi)/2) 

(t-Xi)))  (3.5) 

When  this  signal  at  the  receiver  is  demodulated  with  a  carrier  signal  with  the  same  sweq>  rate 
(e.g.  the  transmitted  signal),  this  results  in  the  signal  conqtonent: 


R«al{^I^Ai**I>(jti')exp(j(-(io^„i„Xi-AtXj^+AXi2/2)) )  (3.6) 


It  is  easily  shown  that  die  frequency  difference  between  a  multipath  conqpoiient  and  die  actual 
transmitted  frequency  dqiends  on  die  delay-time  For  die  i^  signal  die  frequency  diffoence  is 
-Axj,  which  is  always  smaller  dum  the  actual  transmitied  frequency. 

In  case  there  is  tmly  one  padi,  one  dme-dehqwl  vauoa  of  die  transmitied  aignal  is  received. 
Because  of  the  frequent  sweqp,  die  delay-time  cone^ioiids  to  a  certain  frequency-decrease.  Due 
to  rmiltipalh  effects,  signals  widi  different  frequencies  aiiive  at  die  receiver  widi  different 
anqilitude  Aj,  because  die  different  paths  have  different  «nipiitiirtt>  transfer  ftmcripfiy 
Each  frequency  compuueat  of  die  sum  signal  corresponds  with  a  path  detaqt-time.  Tlie  frequency 
differences  can  be  convened  to  delay-time  differences  of  die  differait  ptt««  by: 
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t2-Ti  ”  -  <3. 7) 

Now  from  frie  frequency  q)ectrum  of  die  received  signal  die  nundier  of  padis  (signal 
conqxments),  their  delay-time  and  attenuation  can  be  determined.  In  figure  3.3  an  eunqde  in  case 
2  paths  with  delay-time  difference  are  received,  is  shown. 


Fig.  3.3;  Frequency  chirp  wave  modulation 

Since  no  coherent  detection  can  be  »■>  phase  infosmation  is  avaiUile.  As  for  the 

farowSiand  pulse  measurements,  only  die  energy  of  eadi  frequency  peak  can  be  detennined.  But 
as  distiiict  from  the  broatBauid  pulse  meatmements,  no  qiecial  power  requirenients  for  die 
transmitter  are  necestary.  Although  die  measurement  occiqiies  a  wide  frequency  band,  the 
instantaneous  received  and  demodulated  «igMi  is  a  sum  of  narrow-band  compoaenta  whidi 
energy  can  be  integrated  over  a  sweq^  period.  Hietefine  it  is  not  neceaaacy  to  tranamit  at  very 
high  power  levels  to  achieve  a  reasonable  t^nal-i»-iioise  ratio  at  the  reoeivar. 
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3.3  Pseudo-Noise  sequence  modulation 

This  technique  is  very  close  to  die  broadband  pulse  modulation  technique  as  described  in  3.1. 
However,  two  inqxxtant  problems  are  solved. 

Instead  of  one  pulse,  a  sequence  of  pulses  is  transmitted  on  a  constant  fiequency.  By  coirelatiaa 
techniques,  the  leqionses  of  many  pulses  can  be  summed  into  one  single  pulse.  This  is  possible 
because  of  the  properties  of  the  pseutfaMioise  sequence. 

Bom  the  result  after  correlation  the  delay-times  of  the  different  paths  can  be  detennmed,  similar 
to  die  result  obtained  with  the  broat&and  pulse  measurement  technique.  However,  no  high  power 
transmitter  is  needed  and  by  using  areference  carrier  signal  at  die  transmitter  frequency,  ctdierent 
detection  is  possible  so  that  the  ariqilitude  as  well  as  the  RF-idiase  of  each  path  can  be  measoted. 
Only  the  problem  of  pulse-interference  caused  by  finite  pulse-duration  is  still  left 
Another  advantage  of  this  technique  concerns  the  data  processing.  The  received  signal  is  a 
relatively  long  time-continuous  signal  which  can  easily  be  A/D-converted  and  stored  on  disc. 
Instead  of  real-time  processing,  the  data  can  be  processed  at  later  instants.  This  makes  it  feasible 
to  do  more  complicated  data  processing. 

As  modulating  waveform  for  the  measurement  system,  this  FN-sequence  nxxlnlation  has  been 
chosen,  because  it  has  die  following  advantages  over  odier  waveforms: 

1.  The  transmitted  signal  is  a  low  power  signal. 

2.  The  received  signal  is  easy  to  detect  coherently,  which  makes  it  possible  to  measure  also  the 
RF-phase  differences  between  difierent  paths. 

3.  This  waveform  is  very  useful  tor  digital  processing. 

This  technique  will  be  described  in  mere  detail  in  die  next  diapler. 

It  will  be  diown  there  diat  PN-sequence  modulation  offers  die  principal  possflnli^  to  regain  die 
desired  padi  parameters. 
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4  PSEUDONOISE  SEQUENCE  MODULATION  AS  CHANNEL 

MEASUREMENT  TECHNIQUE 

PN-code  sequeoce  moduUtioD  is  a  land  of  direct  sequence  qxead  qtectrnm  modnlation,  [Dixon, 
1976]. 

In  diis  cluqMer  die  properties  of  a  specific  class  of  FN-sequencea,  Maxinnim  Lengdi  (ML) 
sequences,  are  discussed  and  the  usefulness  of  PN-sequences  as  channel  sounding  wavefoon  is 
shown. 

Further,  the  measurement  system  and  die  design  of  the  transmitter  and  receiver  systems  ate 
described. 

4. 1  Theory  of  Pseudo-Noise  sequences 

The  broadband  channel  measurement  technique  is  based  i^on  the  correlation  properties  of  die 
Pseudo-Noise  (PN)  sequence.  The  sequence  consists  of  code  bits,  also  called  ch^,  wfaith 
assigned  discrete  values  [-1, 1].  The  length  of  a  code  chi^  is  T^  and  Rf. «  l/T^  is  called  die  ch^ 
rale. 

The  length  of  the  sequence  is  K  chqis.  after  which  it  rqieats  itself  again.  The  repetition  time  of  the 
PN-sequence  is  T^  »  KT^. 

There  are  different  kinds  of  PN-code  sequences.  We  have  chosen  for  a  Maximum  Length  code, 
which  has  some  interesting  properties  for  the  use  as  channel  sounding  waveform: 

1.  It  is  easy  to  generate  a  Maximum  Lengdi  (ML)  code  sequence. 

2.  The  auto  correlation  has  a  well  defiiwd  maximom  vahie,  which  is  readied  within  2Tq, 
elsewhere  die  auto  correlation  fimctko  has  a  constant  low  value. 

4.1.1.  A  practical  inqilenientatioa  to  generate  a  PN-aeqnence 

One  property  of  a  Maxhnan  Lengdi  FN-sequence  is  that  conaecutive  bits  of  the  sequence  are 
formed  out  of  the  last  L  tuts.  This  rebtioiiship  holds  even  after  die  last  bit  of  the  sequence  e.g.  die 
*«pK«»ce  starts  over  again.  Tina  9pe  of  FN  sequrnces  can  be  g— with  a  shift  roister 
generalar.  This  generseor  consists  of  a  number  of  shift  legisters  for  slon^  asid  diifthig,  md  a 
number  of  mothilo  2  adders  called  tsys.  Figure  4.1  slimnrs  a  shif>iir|leirrg— «Taiwoftwtgdilf»10 
widi  one  tap  (Hoiiiies,  1982]. 
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Fig.  4.1:  Shift  register  generator  for  generating  an  ML  PN-sequence 

The  output  of  the  modulo  2  adders  is  fed  back  to  die  fost  stage  of  die  shift  register.  Each  clock 
pulse  transfers  the  contents  of  stage  j  to  stage  j-f  1.  The  shift  register  sequence  is  defined  to  be  the 
output  of  the  last  stage.  The  sequence  is  conqiletely  detenmned  by  the  shift  register  leogdi,  the  tq> 
configuration  and  the  initial  contents  of  the  shift  register.  The  actual  contents  of  die  shift  register 
s^.l  a]^.2  ..  a^.L  is  referred  to  as  die  state  of  die  system.  We  can  describe  the  generated  sequence 
by  the  fundamental  linear  recursion  relation,  relafing  the  feedback  sequence  element  a^  and  die 
contents  of  the  register 


•k  ■  «k-l®  =2  •k-2®  ••®  Cl  ajj.L  (4.1) 

where  Cj  e  (0,1 }  denotes  die  feedback  relatxms.  The  symbol  9  stands  for  modulo  2  addition.  The 
vsJues  a^  ai^  j ..  a^  are  termed  the  initial  conditkais. 

For  diese  shift  register  generators,  all  possible  states  except  die  all-zero  stale  occur  during 
sequence  generatioa. 

This  special  group  of  shift  rBgjstw  generators  ve  a  Maximal  Length  generators 

because  the  lengdi  of  the  sequence  equals  the  maximal  nmnfcer  of  different  states  except  the  sll- 
zero  state.  Note  that  foe  all-zero  state  generates  a  of  length  1.  ■«!  is  tmufarwut  after 

each  clock  poise  into  dself  again.  For  a  shift  roister  of  L,  a  maximal  i*«git»  has 

lengdiK  with 


K  -  2^  -  1. 


I 


(4.2) 


I 
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4.1.2  Definition  of  die  auto-cotrelation  fiuction 

To  describe  the  coneladon  properties  of  ML-sequences  we  d^ne  bodi  die  continuons  anto- 
ctmelatioo  function  Rfr)  and  the  discrete  auto-coneladon  function  R(in). 

The  continuous  R(t)  auto-correlation  function  is  given  by: 


R(t) 


(l/KTg)  J  PK(t)PN(t-t)dt 


(4.3) 


The  value  of  this  function  is  [Holmes,  Dixon]: 


1-(K+1) lt|/KT.  for  |t|  <  T_ 

R(T)  -  “  (4.4) 

-1/K  for  |t|  >  Tg 


So  there  is  only  one  triangular  correlation  peak  widi  width  ZT^;  elsewhere  die  correlation  is 
constant  with  value,  -1/K.  A  drawing  of  this  function  is  given  in  figure  4.2. 


Hg.4.2:  Continuous  auto-correlatkai  ftactioo 
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The  discrete  auto-ccneladon  function  R(in)  is  slightly  diffeienL  Let  die  synibols  bj  of  die 
maximal  length  sequence  be  bj  =  -1  or  bj  =  i.  i  =  1^,  ..  JC  and  define  modulo  2  addition 


according  to: 

-1®  -1  - 

1 

(4.5a) 

1®  1  - 

1 

(4.5b) 

-1®  1  - 

-1 

(4.5c) 

1®  -1  - 

-1. 

(4.5d) 

Here  we  use  the  elements  -1  and  1  instead  of  1  and  0  respectively.  Note  diat  because  of  definition 
(4.5),  die  all-zero  state  is  the  state  in  which  all  stages  of  die  shift  register  contains  the  etement  1 
rather  than  -1.  By  definition  (4.3)  the  operation  of  modulo  2  addition  equals  muldplicaticMi.  The 
discrete  auto-correladon  function  of  the  sequence  bj  b2 ..  bj^  is  now  given  by: 

R(m)  -  (1/p)  £bnb{(n+m)  mod  K)' 
n*i 

Since  all,  except  the  all-zero  states  ate  passed,  the  maximal  length  sequence  contains  (2^-l)/2  -  1 
tunes  the  element  1  and  (2^-l)/2  times  die  dement  -1.  For  m  ^  0,  the  auto-correlation  function 
equals  the  sum  of  these  elements  divided  by  the  sequence  lengdi  according  to  (4.2).  For  exact 
overlqi,  m=0,  die  auto-ccnrelation  function  becomes  unity: 

{1  fox  m  -  0 

T  (4-7) 

-1/ (2^-1)  -  -1/K  foe  m  #  0 

The  ouqint  of  die  discrete  auto-ccrrelation  function  is  dnwn  in  figure  4.3. 


Amplitude  Amplitude 
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Pan  of  PN-eeqaence 


Fig.  4.3:  The  disaete  auto-correlation  function 

The  coirelatioa  function  R  has  some  in^rortant  properties: 

The  function  has  one  nuuuiinim  value,  which  has  a  tiiangnlar  shtqre  with  width  2Tg  for  R(t) 
and  a  square  sh^re  with  widdi  T^  for  R(m).  Outside  the  maximum  value,  die  auto- 
coneladon  has  a  low  constant  value  that  dqiends  on  the  wpienry  lengdi.  For  die  FN 
sequence  generator  of  figure  4.1  the  shift  register  has  L  -  10  stages  so  die  ndo  is  2^^-l  = 
1023  >  30.1  dB. 

Takmg  die  comlation  is  a  linear  operation.  Any  version  of  die  FN  sequence,  shifted  in  tnw 
or  mult^lied  in  anqilitiide,  results  in  a  similar  coneladon  fimction  widi  die  «»nm  dififereoces. 
The  sum  of  any  corntrinadon  «rf  N  sequences  will  result  in  a  •  sum  of  correlation 
functioiis. 
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4.1.3  PN-sequence  modulation 

A  PN-sequence,  which  consists  of  code  chips  widi  the  values  [-1,+1],  can  be  BPSK  (Binaiy  Phase 
Shift  Keying)  modulated  <m  an  RF-camer  signal  and  be  transmitted.  This  transmitted  signal  can 
be  writtoias: 


a^(t)  -  PN(t)oos(o>ct)  (4.8) 

The  received  signal  in  a  mult^ath  situation  is  now  found  with  (2.2)  to  be; 


rr(t)  -  Real {^I^Aiexp(j4i’)PN(t-Ti>  exp  ( j  (a>c  (t-Xi) ) ) 


-^X^AiPN{t-Xi)lcos(4i)coBe0ct  +  sin(^i)8ino>gt) 


(4.9) 


where  o)g  [rad/sj  is  the  carrier  fiequency  and  ^  is  die  resulting  RF-pbase  shift  =  ♦£’+  Wctj. 

As  indicated  in  4.1.2,  calculating  the  correlation  of  a  signal  is  a  linear  opeiatioD.  Therefore,  the 
coherently  demodulated  mult^th  signal  can,  after  correlatian,  be  written  as  an  In-phase  signal 
R](t)  and  (}uadrature-phase  signal  RQ(t); 


Ri  (X) 


N 

I  Ai  .R(X-Xi)  .C08(44) 
i-1 


(4.10) 


Fq(X)  -  2  Ai.R(X-^i)  8ln(^i) 
i-1 


(4.11) 


« 


widi  R(x)  is  die  conthmons  conelation  ftmcdon  of  die  FN-aequeace  widi  peak  vabe  atx-  0. 

Now  it  can  be  seen  diat  the  amplitude  and  RF-phase  inftiniiation  is  contamed  in  die  I-  and  Q- 
oulput  signab  Ajcas(bi)  Ai*in(^),  of  the  coberent  demodulator. 

The  above  coasidetatkins  show  diti  conelatioa  of  the  rectived  FN-modnlaied  signal  widi  the 
same  FN-code  offers  die  principal  possibility  to  regain  (he  desired  path  paramelem:  relative 
anqriitnde,  RF-fdiase  and  delay-tiine  of  the  i^  propagation  pndL 


!*■ 

I 


i 
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4.2  Tlie  measurement  ^stem 

In  this  section  the  measurement  set-up  and  die  signal  wavefonns  that  occur,  are  given.  The 
practical  realization  of  the  measurement  system  is  described  in  detail  in  [Fbenstra,  1990]. 

4.2.1  Transmitter 

Figure  4.4  shows  a  block  diagram  of  the  transmitter  which  is  used  to  transmit  die  BPSK 
modulated  PN-waveform.  To  realize  the  transmitted  signal,  we  have  used  a  far  more  conqilicated 
scheme  but  with  diis  figure  it  is  easy  to  explain  die  signal  generation. 


Rg.  4.4:  Block  diagram  of  the  transmitter 


4.2.1. 1  Modulation  of  die  FN-sei]neiice 

The  PN  seipienoe  is  generated  by  die  shift  register  generator  as  given  in  figure  4.1,  with  a  rale  of 
Rq  >  lO^  chips/kecond  .  The  shift  leglslBr  has  L  -  10  stages  so  die  FN  sequence  kagdi  UK* 
2^^-l  >  i023.  I^pon  die  last  bit  die  sequence  siaits  over  agahi.  The  onqmt  of  die  W-seqnence 
generato'  uatitneoaotimioassigaaliadiecdianathiiethscretesignal: 
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PN(t)  -  I  b„  p<t-nTc), 
n-1 


(4.12) 


with  Tf  =  the  symbol  duration  time,  the  symbol  shape  p(t)  is  rectangular  and  given  by: 


p(t)  - 


for  0  S  t  <  Tq 
for  t  <  0,  t  ^  Tj. 


(J.13) 


In  (4.12)  the  element  bg  e  {-1,1 )  stands  fw  symbol  n  of  the  PN  sequence.  The  exact  value  can  be 
determined  with  the  fundamental  linear  recursive  relationship  for  the  PN  sequence  generator  that 
is  used: 


-  bk-3®*’k-10-  (4.14) 

To  Start  the  sequence  the  initial  conditions  bj  b2  ■■  bjQ  must  differ  from  the  all-zero  state.  Note 
that  the  upper  index  of  the  sum  of  (4.12)  is  •»  rather  than  the  length  K  =  1023  of  the  PN  sequence. 
But  since  b^  s  b^  0,0^  K  with  (4.12)  an  infinite  sequence  of  consecudve  PN  sequences  is 
denoted. 

The  time  continuous  PN  sequence  PN(t)  is  BPSK  modulated  on  the  transmitted  signal,  resulting 
in  the  signal  (4.8): 


Bjlt)  -  PN(t)cos(2l(f(.t)  (4.6) 

Fcr’  each  measurement  die  frequency  f|.  of  diis  carrier  is  chosen  in  the  VHF-band.  The  power 
spectrum  of  die  PN-sequence  is  found  to  be  [Ziemer  and  Peterson,  198S,  p.  387]: 


Spi,{f)  -  »(f-Bifg) ,  (4.15) 

widi  Pjj  •  1/K^,  Pg,  ■  [(K+l)/K^]  sii^(niK/K)/(mK/K)^,  K  is  the  sequence  lengdi  and  fg  « 
1/KTf.  Figure  4 J  shows  die  power  spectrum  of  a  FN  sequence. 
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Fig.  4.5;  The  power  spectrum  of  maximum  length  a  PN  sequence 

Because  of  the  periodicity  of  the  sequence  the  spectrum  has  energy  only  at  frequencies  diat  are 
multiples  of  the  reciprocal  of  the  sequence  duration  time  KT,..  The  nuniber  of  fi«quency 
conqtonents  within  one  lobe  equals  the  number  of  symbols  K  widiin  one  period  of  the  PN 
sequence.  Most  of  the  energy  of  the  sequence  falls  within  -Tg*^  <  f  <  Tg'*.  In  our  case  Tg  =  10"®s 
and  thus  most  energy  is  contained  in  a  bandwidth  of  2  MHz. 

4.2. 1.2  Addition  of  a  reference  carrier  component 

To  be  able  to  determine  the  complex  amplitude,  both  in  magnitude  and  RF-phase,  the  signal  has 
to  be  detected  coherently  at  the  receiver.  Therefore  a  reference  carrier  is  needed  at  the  receiver 
which  is  phase  coherent  with  the  transmitted  signal.  At  the  receiver  this  reference  carrier  must  be 
recovered  fimn  the  received  signal. 

Despite  that  the  modulated  PN  sequence  does  not  contain  a  carrier  component,  it  is  still  possible 
to  regain  the  carrier  with  a  Costas-loop  G^emer  and  Peterson,  198S,  p.  26  Iff).  But  in  practice  the 
received  signal  to  noise  ratio  turned  out  to  be  too  low  to  have  successful  carrier  regenentioo  widi 
this  Costas-loop.  Purthermore  tihe  Costas-loop  is  very  sensitive  to  interference. 

For  synchronization  at  die  receiver  we  dierefore  transmit  an  unmodulated  carrier  conqxitient  This 
carrier  is  very  narrow  band  and  can  therefore  be  filtered  out  easily  with  sufficient  signal  to  noise 
ratio  at  the  receiver.  The  filtered  signal  is  then  fed  to  a  {diase  locked  loop  (PLL). 

To  let  bodi  symbols  bj  •»  -1  or  bj  «  +1  have  the  same  energy  upon  recqidon,  die  carrier 
component  is  mad'  in  quadrature  phase  by  a  90  degree  {4iase  shift  iK^th  (4.8)  the  transmitted 
signal  S|^t)  finally  becomes; 
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SjCt)  “  PN(t)cos(oij.t)  +  Acar8in«Bct) , 
with  a>c  =  2itfg  and  ^  is  die  selected  frequency  in  die  VHF-band. 
4.2.2  The  Receiver 

In  figure  4.6  the  block  Hiagnun  of  the  receiver  is  shown. 


Fig.  4.6;  Block  diagram  of  the  receiver 


4.2.2. 1  Demodulatioo 

A  local  carrier  tranqwaes  the  received  signal  r^ft)  to  a  fixed  IF  frequency  band  around  Wjp-  The 
resulting  signal  r,]p(t)  is  filtered  twice. 

A  narrow  band  is  used  to  filter  out  the  added  reference  canfer  congiooent  as  described  in 
section  4.2.I.2.  The  filtered  ovqxit  signal  is  used  by  a  carrier  recovery  networic  (PUL)  to  regain 
the  reference  carrier  signal  This  carrier  makes  it  feasible  to  detect  die  sequence  FM(t)  coherently. 
Because  of  die  multipadi  characteristict  of  die  radio  cfarninel  (be  reference  carrier  that  is  received 
is  the  sum  of  all  carrier  corrqwneats  diat  have  propagated  via  different  padis,  and  dierefare  have 


TNOrapoit 


P»gc 

29 


different  but  constant  RF-phases  and  amplitudes.  So  the  recovered  carrier  has  a  constant  RF- 
idiase  difference  relative  to  die  instantaneous  phase  of  each  padi. 

Beside  the  carrier  conponoit,  the  IF  signal  r^ j[:(t)  contains  die  responses  of  the  FN  signals  of  all 
piiriia  and  a  noise  term.  This  noise  is  assumed  to  be  stationary  white  Gaussian  noise  widi  a  two 
sided  spectral  density  of  No/2. 

We  are  primarily  interested  in  die  determination  of  the  amplitude  Aj,  RF-]diaae  and  delay-time 
tj  of  the  i*  path. 

Since  the  system  is  linear  we  continue  die  system  descrption  therefore  by  describiiig  the 
demodulation  of  one  single  path.  The  total  reponse  is  dien  found  by  addition  of  all  individoal 
path  responses.  The  signal  component  of  path  i  in  the  IF  band  is  given  by: 


-  Ai  Ptr(t-ti)co8(a>ij.t-a>£.Ti-9i')  (4.17) 

-  Ai  PN(t-ti)C08(«)iyt-9i) 

with  (oq;  >>  2i(fjp,  4{'  is  the  RF-jdiase  of  path  i  due  to  reflections  and  odier 

mechanisms  except  padi  delay-time. 

To  limit  the  noise  power,  die  signal  r{p(t)  is  filtered  by  a  2  MHz  wide  filter,  resulting  for  padi  i  in 
the  signal  rgjpift).  Because  of  the  finite  bandwidth  of  diis  filter,  the  spectrum  of  die  PN-sequence 
is  distorted  slightly,  and  also  some  cross-talk  to  die  quadrature  channel  occurs.  The  effect  of 
cross-talk  is  diat  a  signal  udiich  is  in-phase  with  the  reference  carrier,  will  not  only  result  in  a 
signal  in  die  in-iAase  channel,  but  also  in  die  quadtaAite<baniieL 

This  can  be  modelled  in  die  following  way.  When  a  FN-modulaled  signal  in-phase  with  die 
reference  carrier  is  fed  to  die  receiver,  in  die  in-idiase  diannel  die  almost  undistorted  (U)  reidica 
of  die  FN-signal  will  be  found:  FNy(t)  •  FN(t),  whereas  in  the  quadratnte-channel  a  small  but 
heavily  distorted  (D)  FN-signal  results:  FN0(t)  •>  0.  The  exact  descrption  of  die  filtered  signal  of 
(4.17)  widi  cross  talk  conpmients  can  now  also  be  written  as: 


rBijri(t)  -  (PIlB(t-X^)COS(ai2]>t-4^) 

+  P»o(t-ti)ain(tD„.t-4i))  (4.18) 

I 

f. 
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The  signal  rgjpCO  is  fed  to  two  mixen,  driven  by  the  regained  reference  carrier  and  a  90°  (diase 
shifted  version  of  this  carrier.  Since  die  two  carriers  form  an  crthononnal  base  in  signal  qiace,  no 
information  is  lost  by  this  operadiMi.  The  mixers  transpose  die  FN  modulated  signal  from  the  IP 
band  to  die  reqiective  In-phase  and  Quadrature-phase  basdiand  signals  rjj(t)  and  rQjO). 

The  higher  frequency  products  are  filtered  out  by  die  lowpass  filters  widi  inqnilse  response  hj(t) 
and  hqft).  These  filters  ate  almost  equal,  however,  dieir  reqwnses  h]<t)  and  hQ(t)  ate  treated 
differendy  since  slight  differences  in  these  responses  influence  die  data  processing  results.  Now 
let 


Agi  -  Ai/2  coa(4j^). 

(4.19) 

*si  ”  *i/2  ain(4i). 

(4.20) 

PNii(t)  -  PNo(t)  *  hi(t). 

(4.21) 

PNQi(t)  -  PNntt)  *  hi(t). 

(4.22) 

PNlQ(t)  -  PND(t)  •  hQ(t), 

(4.23) 

PNQQ(t)  -  PNa(t)  *  llQ{t), 

(4.24) 

uhere  PN|q  and  PNqj  are  die  cross-talk  conqxments  from  die  I-  to  the  Q-channel  and  vice  versa. 
These  cross-talk  conqxments  ate  very  smaU. 

By  (4.19)  and  (4.20)  the  multqiadi  parameters  Aj  and  f|  for  die  anqilitude  and  RF-phase  of  path  i 
ate  replaced  by  and  A,j,  die  amplitudes  of  die  signal  in  die  in-phase  and  quadrature-phase 
channels.  The  relations  between  Aj,  fj  and  Ag  and  Agj  are; 


*i  ■  +  *8i^>  (4-25) 

♦i  -  arctan(Xgi/Aj.j^)  (4.26) 

Now  the  baseband  ouqiut  signals  are: 

(4.27) 

rQi(t)  -  -A,iPNQQ(t-Ti)  +  AciPHjQ(t-Ti)  .  (4.28) 


The  symbol  *  denotes  die  convohitkm  operation. 
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While  die  cross-talk  conqMnents  are  very  small,  die  essential  terms  of  (4.27)  and  (4.28)  for 
liiidier  cakulatkms  are  those  terms  containing  FNjj  and  PMqq- 

4.2.2.2  Data  san^ling 

The  ou^ts  of  the  filters,  i^t)  and  rgft),  are  fed  to  two  indqiendent  A/D  converters.  Their  ontpnt 
is  stored  on  disc.  In  our  system  we  record  D  *  8  consecutive  PN  sequences  of  K  «  1023  bits  each. 
The  chqi  rate  is  10^  s'^  so  one  PN  sequence  lasts  1.023  ms.  The  duration  of  one  measurement  is 
8.184  ms. 

The  resolution  of  the  A/D  converters  is  8  bils/sanq>le.  This  result  in  2S6  anqdimde  levels. 

For  each  symbol  of  the  PN  sequence  M  ^  8  samples  ate  taken.  Since  the  lengdi  of  die  PN 
sequence  is  K  s  1023,  the  total  number  of  san^les  per  measurement  is  DAf  JC  «  8.8.1023  ° 
65,472. 

We  record  8  PN  sequences  instead  of  one  sequence  in  order  to  be  able  to  interleave  the  8 
corresponding  samples  of  each  sequence.  By  diis  interleaving  operation  die  sample  rate  is 
increased  virtually  by  a  factor  8.  The  interleaving  operation  is  described  in  the  following  chapter. 

If  the  duration  between  2  sanqiles  of  the  A/D  converters  is  denoted  by  T^  die  recorded  sanqiles 
ate  described  by: 

I„  -  for  n  -  1,2, . 

On  "  for  n  -  1,2, .  .,D.M.K, 

Let  us  now  observe  dw  first  2  sanqiles  taken  fiom  PN  sequence  number  1. 

If  die  sanqile  rate  would  exactly  equal  a  mnltqrie  of  the  transmiller  bit  rate  and  in  die  absence  of 
noise,  sanqile  number  I+8K,  die  index  of  die  first  sample  of  die  second  PN  sequence  would  be 
taken  at  dw  same  instant  as  sanqile  number  1 .  The  sanople  rate  R ^  «  T ^  of  dw  converters  is 
chosen  slighdy  lower  fiom  8  times  dw  diqi  rale  R^  ■  1  MHz  (the  exact  sanqile  rate  is  R^  > 
7.999878  MHz).  This  nwfces  diat  equivalent  sanies  of  subsequent  FN-sequences  sre  taken  1/64 
ps  shifted  in  tinw  so  dwt  dw  sanqiles  dw  8  sequences  can  be  interleaved  in  dw  processing 
stage.  This  interieaving  gives  a  virtual  hnease  of  dw  sanqde  rate. 

By  using  diis  techniqne  aw  incieaae  dw  tain|de  rate  viitaally  iqi  to  64 . 1(1^  s'^. 

As  described  beftne,  aw  can  now  add  an  path  leqfwnses.  Togedwr  aridi  dw  noise  tenns,  adiich  are 
indicated  as  n|(nTy^  and  ngfoTy^,  dwse  fyaa  dw  sanies  stated  on  disc: 


(4.29) 

(4.30) 
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H 

+  AaiPNgj (nTjuj-Ti)  +  nj  (nTj„,) ,  (4.31) 
N 

Q„  -  ^I^-AsiPNQQ(n!rM>-'i>  +  +  nQt'iTjuj) ,  (4.32) 


for  n  =  1.2,..JDA1.K.  All  noise  san^les  are  siq>posed  to  be  statistically  independent 


j 


I 


I 


TNOraport 


Page 

33 


3  PROCESSING  OF  THE  RECEIVED  SIGNAL 


In  this  cfaqxer  the  pracetsmg  of  the  received  signal  is  described.  The  signal  has  been  stored 
means  of  die  recorded  samples  on  disc  as  discussed  in  cbqiter  4.  The  processing  is  perfotined  by 
software  on  an  AT  personal  con^wter.  and  consists  of  dsee  parts: 

1.  ralculatimi  of  the  ccerelation  fimcdon; 

2.  Interleaving  of  the  conelatioo  functions; 

3.  Calculation  of  the  multipath  characteristics. 

The  calculations  result  in  die  desired  parameters:  relative  anqilitnde,  RF-phase  and  dday-time  of 
all  relevant  paths.  These  results  may  then  be  transformed  to  the  freipieiicy  domain  to  {aedict  the 
narrow  band  fieqnency  re^onae  of  die  radio  channel,  or  used  to  feed  an  adaptive  equalizer.  When 
a  large  numher  of  channels  are  measured,  also  statistical  infosmation  can  be  gathered  from  die 
data  for  duvacterization  of  different  types  of  terrain. 

This  part  the  data  processing  is  described;  for  mue  detail  on  this  subject  see  ^pendix  B. 

S.l  Calculation  of  die  correlation  fraction 

The  discrete  auto-conelatian  ftuction  of  die  FM  sequence  has  been  calculated  in  dispier  4  and  is 
found  to  be: 


R(«)  -  (1/K)  ,»d  K  - 

f  -l/(2^-l)  -  -1/K  for  B  It  0 
I  1  for  B  ••  0. 


(4.6) 


where  dM  symbols  bj  i-lXX  represent  die  FN  sequence. 

There  are  a  numher  of  reasons  not  to  use  diis  conelatioo  flmctian  (4.6)  straight  forward  in  our 

1.  The  diacrese  conelation  ftmction  (44)  is  a  timfr4isaeie  ftmction.  The  synilKd  halex  m 
conespondswidt  one  drip  duration  time.  In  onrqntMimiiatitw^uT-  iQ-^s. 

The  delqr-time  differences  diat  can  be  expected  are  of  die  same  order  of  majniiwitf  as  the 
widdi  of  the  conelation  ftmdion.  Jf  conelation  ftmction  (4.6)  is  used,  it  is  not  m 

distinpridi  psAs  whh  dday-time  differeiices  wMhin  die  qnhbol  time. 


I 


I 
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Sinoe  we  want  lo  be  abk  to  distingauh  paflia  with  del^r-time  diffcRaces  as  «nan  as 
possible,  8  instead  of  1  (and  after  mtedeaving  64)  conelation  ftmctioa  values  per  symbol  are 
calrulairid. 

2.  The  leqionse  of  an  occnmng  propagatiaa  padi  is  contained  in  die  samples  from  die  I-fbase 

as  well  as  the  <}-phase  channeL  The  !•  and  (^-san^les  are  processed  separately.  They  contain 
die  RF-fdiaae  infonnatimi. 

The  duratiao  of  the  coerelatioa  peak  of  each  padi  should  be  as  short  at  possible,  in  order  m  have 
minimum  overlap  between  die  correlation  maxima  belonging  to  different  padis  widi  almost  equal 
delays.  The  length  of  this  reaponse  partly  dqiendt  on  die  fimction  that  is  used  to  correlate  the 
received  samples.  One  of  the  nqiut  ftmctions  tor  the  correlation  operation  is  the  sequence  of  the 
received  sanqiles.  The  odier  input  fimction  can  be  dioaen  so  diat  die  duration  of  the  ouqmt 
correlation  maximum  is  as  short  as  possible. 

The  fimction  diat  is  used  heie,  is  dwsen  slightly  different  from  the  FN-code  signal  (which  would 
result  in  a  trisngnlar  conelmkm  peak  diat  is  wide).  The  signal  has  die  same  value  asdieFN- 
signal  (-fl  or -1)  for  the  first  sample  period  of  a  ch^  (T  >  1/8  T^  and  is  tero  elsewhere. 

This  fimctirm  can  be  deacribed  by; 


“  *>  r(k)/MDl  ’'((k-X)  mod  MD)  ^  (5.1) 


widi  V|(  the  discrete  Dirac-fimction: 


1  for  k  -  0 
0  fox  k  #  0. 


(5.2) 


The  fimction  Txl  denotes  dm  smaOeat  integer  greater  Oan  or  equal  to  X  and  b^  a  (-1,1)  is  die  PN 
sequence  of  (4.14). 

The  wxUi  of  the  conelation  peak  in  this  case  is  1  IIS.  This  is  half  the  widdi  of  the  ccnelatian  peifr 
of  the  analog  condatian  ftmctian,  where  the  width  of  the  peak  is  2  IIS. 

Both  ftmctiaiis  snd  die  condition  fimction  see  depicted  in  figure  3.1. 


i 

I 


Conelitiaa  bmctkn  wilh  ad^Med  wavefocm 


Rg.  5.1:  Signal*  and  coirdatioo  imictioD 

The  calculation  of  the  ccaelation  functions  of  die  received  !•  and  Q-  signals  can  now  be  described 
fonnally  as: 


Rjln]  -  ^_^I[n+kl  x(k], 
RQln]  -  ^  I  Q{n+kl  x{k]. 


(5.3) 

(5.4) 


wifli  die  nootiaa  z(kH^  uaed  to  denote  IJk].  QM.  and  xM  defined  in  (4.31),  (4.32)  and  (3.1) 
wsnectlvely. 
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To  clarify  the  correlation  an  exan^le  for  a  amplified  case  is  calculated.  In  the  Mtwpliff'^  case 
only  one  padi  is  received  and  die  received  samples  exactly  equal  die  transitted  FN  sequence.  If  we 
assume  that  the  received  signal  and  die  regained  carrier  are  in-idiase,  we  nu^  ftndier  disregard  die 
quadrature  channeL  The  received  sanqiles  sinqily  become: 

I[k]  -  b[  fk/lffil  ].  (5.5) 

Widi  (5.1),  (5.2)  and  (5.3)  the  correlation  function  is  calculated: 


Rjtnl  -  I  x(k)  I[n+kJ 


k— « 


^  Z  birk/Mol]  V[(k-1)  mod  MD]  b  [F  (k+n) /MdI]  « 


-  ^  Z  blj+1)  blPj+d+nl/MDl] 


Z  b[j+l)  b(j+l)  for  0  S  n  S  MD-1 
J— ~ 

m 

Z  bIj+11  blj-lf  (l+n)/MDl)  for  n  <  0,  n  >  MD-1 

j— 


(1  for  0  S  n  S  ND-1 

-l/(2^-l)  for  n  <  0,  n  >  MD-1. 


(5.6) 


The  conrelation  fuDctioo  is  an  nact  re|dica  of  one  symbol  leqiooae  of  the  sequence  I[k]. 

Thia  conespoodence  is  explamed  as  follows:  The  received  conaiat  of  consecutive  symbol 

reaponaes  of  die  FN  sequence-  By  hddng  die  cosielatMai,  samples  ate  added  and  sobnacied.  In 
case  die  received  amnfdea  and  die  odier  hqkit  cooelation  function  overiap,  values  all 
contribute  with  die  same  sign  e.g.  positive  signed  — are  wtHfi  hk)  negative  mpoxt 
mt  sabcraded.  The  result  is  an  ani(dified  value  of  one  symbol  response  vabe. 

In  reality,  the  received  sanqiles  win  not  equal  die  PN  sequence  at  in  (5.5)  but  diey  win  be  spread 
in  time  herunae  of  handwiddi  limit atiuns.  bi  die  resnltiiig  comlatian  function  a 
apnad  is  found. 
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S.2  Interleaving  of  tiie  conelalion  functions 

The  I-  and  Q-  signals  are  sampled  over  a  period  of  8  consecutive  PN-sequeaces.  Howeva,  die 
moments  of  equivalent  sanqiles  in  two  successive  sequences  are  deliberately  chosen  not  to 
coincide,  as  described  in  cluqiter  4.  They  are  shifted  over  1/8  of  the  sample  period  (1/64  ps).  The 
same  is  true  for  the  coitelation  fbnctions.  So  the  8  consecutive  correlation  functions  ate  all 
slightly  shifted  over  1/64  ps  within  one  sample  period.  Now  the  ccmelation  functions  of  die  8 
sequences  belonging  to  die  same  channel  can  be  interleaved.  This  operation  is  shown  in  figure 
5.2. 


Sample  moments  in  a  chip  duration 


Fig.  5.2:  Sanyling-  and  interlesving  operation 

Sj3  indicates  die  fifth  sanqde  of  die  cotreladoo  fimetion  (die  duratka  <rf  srtndi  is  about  Tg)  of 
PN  sequence  1,  $23  is  the  equivalent  sample  of  PN  sequence  2, ,  etc. 

This  results  in  an  equivalertcoirelmion  fimetion  wife  an  8 
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To  remove  die  noise,  tbe  resulting  conelation  function  is  filtered  widi  a  moving  average  filter 
-with  length  8.  The  processing  gain  Pq  ('arhich  is  die  inqnovemeat  of  die  signal-to-noise  ratio)  of 
the  coneladon-  and  interleaving  operadons  together  is  101og(8.1023)  =  39  dB. 

Of  course  this  can  only  be  perfonned  if  the  channel  parameten  don't  vary  in  die  observed  time 
intervaL 

The  advantages  of  this  operation  are: 

A  higher  virtual  sample  rate,  resulting  in  an  8  times  better  time  resolntioii. 

The  system  error  is  now  1/128  ps  =  -f/-  7.8  ns. 

The  signal  independent  noise  contained  in  tbe  sanqiles  of  different  sequences  is  uncorrelated 
and  can  be  filtered  out 

The  new  correlation  fiinction,  after  interleaving  and  filtering,  is  used  to  calculate  die  desired 
multqMth  parameters. 

S.3  Calculation  of  the  multqiadi  characteristics 

Calculation  of  the  multipath  characteristics  firom  die  I-  and  Q-conelation  signals  is  quite 
complicated,  and  therefore,  is  given  in  full  detail  in  ^ipeDdix  B.  Here,  only  a  qualitative 
description  of  these  calculations  is  given. 

The  iiqiuts  for  these  calculations  are  the  I-  and  Q-coirelation  functions  and  the  reference  I-  and  Q- 
corielatiott  functions  (Refj  and  Re^  Refj  contains  information  fiom  PNq  and  PNqj,  Re^ 
contains  infonnation  from  PMqq  and  FNjq).  The  reference  conelatioa  fimctioos  are  determined 
by  the  transmitter  and  receiver  parts  of  die  system  only.  They  are  not  dqieiident  on  die  radio 
channel  c.q.  the  muUqiath  characteristics.  These  reference  coiielatioo  functions  which  coatain  die 
influence  of  die  transmitter  and  receiver  characteristics  due  to  filtering,  phase  shifting  etc.,  ate 
used  to  estimate  die  nailtqiadi  characteristics.  They  are  measured  laboratory  conditioiis 
with  transmitter  and  receiver  ditectfy  connected  togedier. 

The  calculated  coccelation  fonctions  of  die  received  I-  and  Q-tignala,  RjCn]  ud  RQ[n],  are 
detenninedby; 

die  mnhqNlh  diaracteristics  N,Ag,A^  and  1}  for  i  -  1,2,.. 
die  noise  terms, 

as  given  in  (S  J)  and  (3.4)  with  (4.31)  and  (4.32).  Ij  is  die  discrete  delay-time  valiie. 
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‘  By  chosing  die  right  anqilitudes  and  A^j,  and  delay-time  1^  for  path  i,  and  assign  these  to  Refj 

and  Re^,  the  newly  fmned  Refj(i)  and  RefQ(i)  can  be  matdied  witii  the  received  I-  and  Q- 
coirelation  signals  belonging  to  patii  i. 

We  can  now  define  the  error  criterium  of  the  data  processing  to  estimate  the  multipath 
characteristics:  determine  the  number  of  patiis  N  and  for  each  path  its  amplitudes  A^,  Agj  and 
delay-time  Ij  so  that  the  energy  contained  in  the  remaining  sanqile  values  in  tiie  I-  and  Q-channels 
is  minimum.  In  other  words:  find  N,  A^,  A^j  and  Ij  fi>r  i>l,2,..J>I  so  that 

a2  -  X  {  Rjln]  -  X  AjBefitn-li]  + 
n«-««  i-1 

•»  N  _ 

‘  +  I  {  Rflln]  -  X  AQRefjjtn-lil  (5.7) 

n»— oo  i“l 

t 

is  minimum.  Since  tiie  energy  of  Rj[n]  and  RQ[n]  is  concentrated  in  die  peaks  of  the  function 
only,  the  summation  of  (S.7)  can  be  restricted  to  a  limited  number  of  sanqiles.  This  number  must 
be  determined  first  and  will  differ  for  each  radio  channel. 

An  iterating  algorithm  is  used  to  determine  the  multipath  characteristics.  The  block  structure  of 
the  algorithm  is  shown  by  the  flowchart  given  in  figure  5.3. 
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Fig.  5.3:  Flowchart  for  the  rletermiiiation  of  the  mitnber  of  paA  aixl  the  path  parameters 

The  algorithm  starts  with  Nsl  path  arid  r^tthnizes  the  amplitudes  A^]  and  Agj  and  the  delay-time 
Ij.  With  these  parametm,  which  ate  assigned  to  the  reference  conelation  functions  Refj[l]  and 
Refqll],  die  error  I?  of  (5.7)  is  calculated.  This  error  function  is  used  to  optimize  die  assigned 
path  parameters. 

When  die  parameters  are  opUmnm  and  die  maxiimim  value  of  die  enur  functioa  A^,  still 

exceeds  a  chosen  threshold  level,  the  energy  of  the  error  function  is  sqipoaed  to 

originate  from  anodier  padi.  In  diis  case  another  path  is  added  and  the  algorithm  starts  over  again 
to  find  the  optimum  sohdion  for  bodi  padis  together. 

Padis  are  added  until  die  resulting  maximnm  error  beemnes  smaller  than  The 

delay-dme  of  die  new  extra  padi  is  determined  by  die  maxhnnm  of  die  error  function  A%n]. 

The  newly  added  padi  causes  a  decrease  of  the  error  froctian  A^  if  die  energy  in  originates 
from  other  paths  of  die  radio  channeL  In  case  the  energy  of  the  error  ftmedon  is  caused  by 
reemved  noise,  the  dwnge  of  the  error  I?  is  undetermined,  hi  pracdoe  it  turns  out  diat  dw  error  e? 
always  decreases.  Hence  we  may  state  that  die  algorithm  converges. 
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The  last  step  to  obtain  all  mult^ath  characteristics  is  die  detennination  of  the  anqilitude  and  RF- 
pbase  differences  of  the  found  paths  from  the  anqilitudes  and  Agj.  These  are  already  given  in 
(4^5)  and  (4.26): 


I 


+  Aai^) 

(4.25) 

«  arctan(Aai/Aci) . 

(4.26) 

i  ■ 
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6  ACCURACY  OF  THE  MEASURED  PARAMETERS 

Rom  die  channel-measarenient  results  die  relative  amplitude,  RF-phase  and  delay-time  of  die 
different  mult^iath  rrignais  are  detennined  by  processing  of  the  received  signal  as  is  presented  in 
the  previous  chqiters.  How  precise  die  characteristics  of  the  radio  channel  can  be  estimated, 
depends  on  die  accuracy  of  the  measured  padi  parameters. 

In  this  chqiter  die  accuracy  of  the  measured  parameters  will  be  evaluated.  The  sources  of  enor  for 
the  parameters  are: 
an^litude:  noise 

delay-time:  noise,  sample  frequency 

RF-phase;  noise,  phase-fitter  of  the  reference  carrier. 

In  the  following  at  first  the  signal-to-noise  ratio,  and  the  influence  of  phase-jitter  of  the  reference 
carrier  will  be  discussed.  Then  the  effects  of  the  error-sources  on  the  accuracy  of  die  path 
parameters  will  be  evaluated. 

6.1  Signal-to-Noise  Ratio 

The  total  received  signal  crmtains  a  signal  part  S,  and  a  noise  part,  (r  •>  received).  The  SNR 
(Signal-to-Noise  Ratio)  of  the  received  signal  is  defined  as  SNR  S/N^. 

The  SNR  of  the  received  signal  after  A/D-conversion  is  determined  by: 
the  SNR  of  the  received  signal  S/N^, 

the  quantization  noise  caused  by  A/D-conversion  Nq  (q  <=  quantization). 

-  Influence  of  quantization  noise 

The  influence  of  Nq  can  be  determined  by  conqiating  the  SNR  before  and  after  A/D-conversion. 
A/D-conversioa  is  performed  by  an  8  bits  converter.  This  means  diat  die  voltage  range  is 
subdivided  into  2S6  steps.  When  die  voltage  range  is  used  optimum,  and  N,.  is  assumed  to  be  0, 
die  maximum  SNR  that  can  be  reached  is,  [Cariscm,  1975]: 

SW*niax  ■  “  10*lo» (3-25«2)  -  53  dB.  (6.1) 

Because  N^  and  Nq  are  independent,  die  total  noise  power  after  A/D-conversion  is  N|o|  > 

Nf  +  Nq. 


Now  it  caa  be  shown  easily  Aat  die  influence  of  Nq  on  die  total  SNR  becomes  significant  onfy 
when  the  SNR  of  die  received  signal  is  about  SO  dB  or  higher,  hi  practical  field 
measuremrats  such  high  S/N^  will  not  be  common  practice.  More  reasonable  values  are  S^^  <  20 
dB.  Then  the  influence  of  Nq  is  neglectable,  and  die  SNR  after  die  A/D<aaverter  is 
qjproximately  equal  to  S/N^. 

-  SNR  after  signal  correladcm. 

Up  to  now  only  the  SNR  before  signal  processing  (correlation)  has  been  investigated,  however, 
the  SNR  will  increase  significantly  after  correlation.  The  effect  of  correlatiao  and  interleaving  of 
8  PN  series  results  in  a  processing  gain  Pq  of  Pq  =  101og(8.1023)  ■>  39  dB.  So  the  SNR^  (SNR 
after  correlation)  of  die  signal  which  is  used  to  determine  the  pulse  position,  is  related  to  S/h^  by: 

SNRc  -  S/THj.  +  39  dB.  (6.2) 

In  this  result  the  influence  of  die  processing  gain  is  clearly  demonstrated.  Even  with  a  low  S/N^ 
the  SNR^  can  be  quite  good. 

-  Determination  of  the  SNR^;  for  a  single  multqiadi  conqxment 

For  a  given  pulse  shtqie,  the  dme  accuracy  is  determined  by  die  SNR^  after  correlation  as  will  be 
shown  in  paiagnqih  6.3.  SNR^  has  to  be  determined  for  a  ^gle  pulse  (or  path).  For  a  single 
pulse  S  and  N  can  be  calculated  as  follows: 

Tp 

s  -  (Tp)-l  /  |P(t)|2  dt  (6.3) 


1  f“  5 

N  -  (T„)  1  n*(t)  dt  (6.4) 


IfcreTp  is  die  bme  interval  of  a  die  pulse,  Tjj  is  a  time  interval  widi  no  signal  energy  preaenL 
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6.2  Errors  caused  by  phase-jitter  of  the  refereace  carrier 

In  the  signal  processing  phase,  8  successively  received  FN  sequences  are  cmelated 
independently.  After  correlation  the  results  are  interleaved.  In  the  ideal  sibiation  a  smooth  curve 
wiU  result  widi  sample  values  at  1/64  ps  distance.  However,  in  practice  it  points  out  that  die  phase 
of  the  reference  carrier  is  not  stable.  It  contains  a  significant  amount  of  jitter.  This  causes  diat  dw 
amplitude  of  the  signals  in  die  in-phase  and  quadrature-phase  channels  of  die  receiver  are  not 
stable  in  time. 

An  exanqile  of  the  resulting  correlation  peak  afto-  interleaving  is  given  in  figure  6.1. 


} 


ng.6.1:  Cone]at^pmdc,a8.caicnlated  from  measurement  tesiilli,widi  die  eSbct  of  phase-jitter 
(I<  Q)  and  alter  moving-average  filieriiig  Q,  Q) 

By  using  a  moving-average  filter  widi  kngdi  8,  the  variations  are  filtererl  out  (figure 

6.1).  By  sunulalion  die  efitects  of  these  operations  finteriemrmg  with  ihssf  jllii  i  and  filleij^) 
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have  been  evaluated  for  diffezeat  peak-^teak  phase-jitter  values.  It  pointed  out  ttiat  phase-jitter  of 
the  reference  catrier  causes  an  attitude  errin',  but  no  RF-idiase  error. 

6.3  Delay-titne  etrar 

6.3.1  Delay-time  emns  due  to  noise 

In  the  case  of  a  radar  signal,  the  time  t]apscd  between  the  moment  of  transmission  of  the  radar 
pulse  and  reception  of  the  leflected  signal  determines  die  range  to  die  origin  of  die  refkctkai.  In 
general  the  received  signal  is  a  dishMted  replica  of  die  transmitted  pulse  eonbedded  in  noise.  This 
causes  that  the  measurement  of  the  delay-time  has  not  infinite  accuracy. 

The  error  in  die  delr^-dme  of  die  cmrelation  peak  enibedded  in  noise  in  our  case  is  equivalait 
with  the  error  in  the  echo-time  of  a  radar  signal  in  noise. 

Therefore,  calculations  originatiiig  fiom  the  radar  science  [Woodward,  19S3;  Barton,  1969]  can 
be  used  to  determine  die  accuracy  of  the  measured  relative  delay-times  between  the  correlatiaa 
peaks  of  the  multqiath  con^onents  and  die  direct  signal. 

The  delay-time  accuracy  that  can  be  achieved  depmds  on  the  SNR  and  die  shqie  of  the  pulse. 
The  accuracy  will  increase  with  increasing  SNR  and  widi  inaeasing  rise-  and  fidl-times  (whidi  is 
equivalent  with  a  larger  spectral  bandwidth  of  die  signal)  of  die  pulse. 

A  measure  for  the  accuracy  is  die  RMS  (Root  hfean  Square)  error  [s].  For  a  given  pulse  slufie 
the  minimum  RMS-etror  is  given  by  o^,  [Woodward,  1953;  Barton,  1969]: 

0.5  -  l/(PVsNF„)  .  (6.5) 

Ihe  received  pulse  energy  is  maximized  when  the  receiver  characteristics  are  matdied  to  the 
pulse  (matched  fillet  receiver).  For  a  matched  filter  receivct  SNRg,  >  2E^0  ^  ^  ^ 

is  the  energy  in  die  received  pulse  and  Nq  die  qwcttal  noise  denai^.  The  variable  P  contains 
information  about  die  pulse  shape,  it  eqiecuDy  dqiends  on  the  rise-  and  fall  times  erf  die  pulse.  P 
Ugivettby: 
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/  I  3/dt  P(t)  |2  dt 

— OO 

p  -  - 

r I  p(t)  |2  dt 


(6.6) 


P  is  the  RMS  baoidwidth  of  the  pulse  signal,  P(t)  is  the  complex  pulse  shape. 

To  be  able  to  calculate  die  delay-tiine  enor  due  to  noise,  die  SNR  of  the  pulse  signal  has  to  be 
known.  The  receiver  which  has  been  developed  to  perform  practical  measurements,  is  an 
qiproximadon  of  the  matched  filter  receiver.  Therefore  the  received  SNR^  is  sub-optimum,  SNR, 
<  SNR^.  hi  the  following  it  will  be  assumed  that  SNR,  >  SNR^^  to  compensate  for 
mismatdiing. 

In  figure  6.2  the  calculated  standard  deviation  o,  of  the  delay-time  error  is  given  as  function  of 
SNR^  (SNR  a&er  correlation  as  discussed  in  paragtsph  6.1)  for  the  same  pulse-shqie  as  is  used 
in  the  measmement  system,  for  which  P  =  1.86'10^. 


0  M  «  so 

■NR  «  smHMS  Mem  m 

Fig.  6.2;  The  calculated  standard  deviation  of  the  delay-time  error  as  ftaiclioa  of  d>e  faipnt-SNR 
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6.3.2  Inftoence  of  die  meuarement  method  on  die  del^-liiiieencc 

The  effective  lanqile  >»«»«  of  the  ngnal  niker  cooeladon  and  interkeving  is  1/64  ns  >  1S.6  ns.  The 
inherent  enor  caused  by  diis  disaetizatioD  is  +/-  7.8  ns.  Because  the  distiibatiaa  of  the 
discretization  error  is  uniform,  die  standard  deviation  due  to  this  error  is  4.S  ns. 

From  fig.  6.1  it  can  be  seen  diat  for  SNK^  >  45  dB  (iiqiot-SNR  >  6  dB),  the  standard  deviation  of 
the  delay-time  error  is  smaller  than  die  system  standard  deviation  of  4.5  ns.  Bodi  dieae  error  are 
indqiendent,  so  die  total  variance  is  found  by  adding  both  variances. 

For  the  delay-time  between  two  cotrelatioo  peaks  (1  and  2),  the  variance  in  die  relative  delay-time 
is  determined  by  the  variances  in  die  delay-times  of  bodi  peaks  (the  errors  ate  indqiendent): 


6.4  Anqilitude  error 

In  the  chosen  measurement  set-np  there  are  two  causes  of  errors  in  die  measured  amplitude  of  a 
radio  padi. 

1.  An^iitude  errors  due  to  noise.  The  presence  of  die  noise  will  cause  a  statistical  deviation 
that  dqiends  on  die  SNR. 

2.  In  die  nwasurement  set-up  the  r^etcoce  carrier  is  recovered  by  means  of  a  Phase  Locked 
Loop  (PLL).  The  refavnce  carrier  fimn  the  PLL  still  has  some  jJiaic  jillrr  Dqieadhig  on 
the  amount  of  jitter,  dds  also  causes  amplitnde  errors. 

6.4.1  Amplitude  error  due  to  noise 

A  carrier  signal  |dns  noise  utedi  is  bandpass  filtered  can  be  written  as  a  sum  of  an  ii>pliase  and  a 
quadrature  phase  oompoocnt 


a(t)  -  (A  +  n^tt)  )  cos(set)  -  nq(t)  sin(ttt)  - 
-  *(t)  cos  (set  ♦  6(t)) 


(C.S) 


In  diis  fatmnla  R(t)  and  t(t)  arc; 
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R(t)  -  [{A  +  ni(t))2  +  (nq(t))2)0-5  (6.9) 

4(t)  -  arctui((nq(t)  )/(A  +  n^(t)))  (6.10) 

Thu  equatkn  relates  tiie  envelope  and  phase  of  a  earner  bandpass  ww  to  die  signal  and 
noise  components.  The  probability  density  ihnetion  (PDF)  for  R(t)  and  4(t)  is  given  (Carison, 
1975): 


H  -  . 

P(R,4)  -  -  •3q>(-(R2  -  2ARCOS  4  +  a2)/2H)  (6.11) 

210) 

This  is  a  joint  PDF,  and  R  and  4  are  statistically  dependent. 

For  large  SNR  >  A^/2N  (ahere  N  is  die  noise  power  in  the  pass-band  of  the  filter),  die  probability 
density  function  for  die  amplitude  can  be  estimated  by: 

Pr(R)  -  (R/2)CAK)1/2  ,xp{-(R-A)2/2H)  (6.12) 

Now  it  follows  for  the  mean  deviation  and  variance  of  the  amplitiide  (Chrisen.  1973): 


R  -  A  (6.13) 

Or^  -  K  (6.14) 

With  (6.14)  die  standard  deviation  of  die  mpHtinfa  error  be  calcnlated  as  fhoctian  of  SNR^ 

In  figure  6.3  die  vahie  (A  +  og)/A  » 1  + 1/  ift'SNR^  u  given  as  fimetion  of  SNR^. 


SNR  S(  oimlMM  HgM  MB) 

Hg.  6.3:  The  calculated  staixUid  deviJttUm  of  the  amiditiide  enor  due  to  noise 

6.4.2  Anqtlitade  enor  due  to  jduue-jitter  of  die  lefeRace  earner 

As  discassed  in  paiagrqih  6.2,  the  pbMte-fiOa  of  the  leAreaoe  earner  mfinMinM  liie  signal  that  is 
received  of  die  8  successively  transmitted  FN-seqnenoea.  anmlation  showed  »♦««*  Hii« 

causes  an  aniplitutle  enor  that  depends  on  the  peak*peak  phase-jilter.  The  calcnlaled  amplittide 
error  is  given  in  figure  6.4.  In  practice  the  peak-peak  phase-jitter  (rf  dw  lefiaence  catrier  will  be 
within  45  degrees,  than  the  maximnin  enor  is  within  0.S  dB. 


ptik-PMk  phaM  IMr  (dagr) 


Fig.  6.4;  The  effect  of  phase-jitter  on  die  amplitode  enor 

6.S  RF-phase  error 

The  calculated  RF-pfaaae  related  to  die  signal  cootent  of  die  I-and  Q-cfaanneb,  will  deviate  when 
noise  is  present.  From  die  ghren  joint  PDF  of  die  ampHhide  and  phase  as  given  in  paragrqpta  6.2.1 
,  the  PDF  (Robatnlity  Density  Function)  of  die  phase  errer  can  be  estiinated  for  SNR  »  1  by; 

p^(t)  -  (a2/2iW)1/2  cos^  »^n^) /2V}  (6.15) 

From  diis  relation  it  foDows  for  die  mean  enor  and  die  variance  diat: 


TNOraporl 


Page 

51 


In  the  way  as  ghrm  in  paragn^h  6.4.1  die  standard  deviation  of  the  RF-phase  is  calcnlated 
as  function  of  SNR^.  Hie  results  are  given  in  figure  6.5. 

50 
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SNB  o<  oofTtimd  tQna  (dB) 

Hg.  6^;  The  calculated  standazd  deviation  of  die  RF-{dM8e  error  doe  to  noise 

For  SNR^  >  20  dB,  die  standard  deviadm  of  die  jdiase  enor  <  S  degr. 

6.6  Conclusion 

The  errors  due  to  noise  are  vety  For  a  SNR^  >  23  dB  (S/I^  -  -14  dB),  die  standard 
deviatkMt  are: 

tune  :22  ns 
amplitade  :  OJdB 
RF-phaae  :  3  degr. 

In  general  die  a^mt-SNR  win  be  much  higher  so  diat  the  enots  due  to  noise  become  very  smaU. 
With  reflect  to  die  calculation  of  die  dehy-time,  the  systematic  error  is  widiin-ty- 7.8  ns. 

Phase  jitter  in  the  reference  canier  canses  only  an  amplitiKle  enor  and  no  phase  error.  Thk  error 
is  not  dependent  on  the  RFpbaae,  h  only  dqiends  on  die  pedk-peak  phase-jitter.  For  peak-peak 
phaae-jiiter  ■  30  degr.  the  ampliRide  error  is  0.2  dB. 
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7  FREQUENCY  RESPONSE  OF  THE  RADIO  CHANNEL 

In  this  cluq>ter  the  method  that  has  been  used  to  calculate  die  frequency  reqxnue  of  a  multipath 
radio  channel,  will  be  discussed.  The  parameters  of  the  single  paths  in  die  radio  channeL  udiich 
can  be  measured  and  calculated  as  described  in  die  previous  duqiters,  are  input  for  the 
calculations. 

One  of  the  problems  that  will  be  addressed  is  die  influence  of  die  unceitainty  in  die  measured 
delay-times. 

7. 1  The  frequency  reqionse 

In  chapter  2  the  model  of  die  muldpadi  radio  channel  is  discussed  in  detail  There  the  iirqmlse 
lesponse  h(t)  of  the  radio  channel  (2.1)  is  given: 

h(t)  -  ^S^hi<t)  -  (2.1) 

where  Aj  is  the  amplitude,  '  the  RF-phase  and  tj  the  delay-time  of  path  i.  N  is  die  number  of 
paths.  The  frequency  response  H(f)  can  be  calculated  by  talcing  the  Fourier  tiansfonn  from  die 
inqiulse  response  h(t).  H(f)  is  given  by: 


H(£)  -  rh(t)exp(-2Kjft)dt 

-•O 

N 

"i-i  • ) ) 


(7.1) 


H(f)  is  a  complex  function  and  is  tfaerefose  described  by  its  anqilitade  and  phaae.  In  stead  of  die 
arrylitride,  we  will  use  here  die  power  lyerlrum  S(t)  >  IH(f)t^.  Widi  (7.2)  the  power  qwetrum  S(f) 
can  be  calculated  from  die  padi  paraurters. 
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S<f)  -  (  i  Alcoa <2ie£Ti-^i*) 


+  {  I  Ai8ln<2*fti-^i’)  )■“ 
1-1 


(7.2) 


The  phase  response  as  fhnction  (rf  frequency  is  given  by: 


N 

4<£)  -  »rct*n({  L  AiC08<2KfTi-4i’)  )/ 

N 

{^I^Ai8ln(2*fTi-4i’>))  (7.3) 

From  (7.2)  we  can  see  that  the  power  q>ectnim  is  firequency  dependent  because  of  the  multipath 
effect  The  variation  in  the  received  signal  power  caused  by  interference  is  determined  by  die 
anqilitude  A|  and  delay-time  of  each  path.  The  {diase  reqxinse  ^f),  as  given  in  (7.3),  also 
deviates  from  the  ideal  linear  case  due  to  multipath. 

For  frequency  management  purposes  die  power  spectrum  S(f)  contains  all  necessary  infcimation 
for  narrow  band  signals.  For  this  type  of  signals  phase  variation  over  the  frequency  band  occiqned 
by  the  signal,  will  be  small.  Therrfote  in  die  following  only  the  power  qiectrum  S(0  will  be  taken 
intoaccosmt 

7.2  Influence  of  the  delay-time  accuracy 

When  calcalating  the  frequency  teqioose  from  the  path  paramwnt  die  {dnue  at  at  certain 
frequency  is  detetnuned  by  die  measured  RF-pbase  at  die  measurement  frequency  and  die 
delay-time  ij.  The  accuracy  (rf  is  of  great  importance  for  the  accuracy  of  S(f). 

If  we  consider  path  J  widi  delay-time  Xj  the  RFfhase  of  this  padi  is  periodical  widi  lAtj.  The 
contribution  of  dna  padi  widi  Aj  in  (7.2)  win  changed  from  its  mhiiiiinm  (-Aj)  into  its 

maodmum  value  (-f  Aj)  by  a  diange  of  die  phase  by  n  rad. 

The  uncerlaiiiQf  of  Ae  delay-times  has  two  maia  causes: 

1.  The  Ante  sampling  rate  ofR,- 64*10^  rhesuses  a  discietizationenorofAt  re  7.8125  ns, 

2.  Noise. 

At  die  frequency  f  re  SO  MHz  a  phase  diift  of  a  rad  is  cansed  by  a  dehQr-time  shift  of  re  10  na. 
This  means  diat  the  uncertainty  in  die  meaanred  delay-times,  whidi  is  in  die  order  of  20  ns,  is 
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about  one  period  of  tbe  RF-signal  in  the  frequency  band  30  -  70  MHz,  ^^lae  the  measmenients 
are  earned  out 

In  the  foUowing  paragraph  the  effects  of  the  delay-time  inaccuracies  will  be  shown  and  a  method 
will  be  psesented  to  partly  con^rensate  for  Oiese  inaccuracies  by  changing  the  delay-times  tj  of  the 
preqtagation  paths  slightly  so  that  the  total  RF-phase  shift  at  the  measurement  freqimicy  is  exactly 
the  measured  RF-phase 

7.3  Calculatioa  of  the  frequency  response  from  the  measured  multqnth  parameter 

It  is  assumed  that  die  measured  delay-time  of  path  i  has  a  Nmnal  distribution  around  the  actual 
value  with  mean  value  Xj  and  variance  from  the  radar  analogy,  as  described  in  chapter  6 
[Barton,  1969],  diis  is  reasonable.  Another  assumption  that  is  made,  is  that  die  variance  of  die 
delay-tiine  is  constant  for  all  paths,  a-^  » <>2. 

The  delay-time  distributioD  must  be  taken  into  account  for  the  calculation  of  the  power  plectrum. 
The  distribution  function  for  die  delay-time  of  a  single  path  is  given  by: 


•  l/VtaWT^l  •xp{-(t-Ti)2/'(2o2)  )  (7.4) 

Each  delay-time  Tj  is  a  stochastic  variable  now.  Widi  (7.4)  we  can  ralmlafe  die  probability 
density  function  p((t)  of  each  ttelay-time  set  (T],‘^,..,t^)  chosen  in  die  N-dimensional  time  space: 


N 

“  Ptitr  • ’T*  “  n  Pxi  (t)  (7.5) 

Because  of  the  uncertahity  in  die  delqr-tmies  of  die  padis,  also  die  power  tpectmm  at  a  tvwfin 
frequency  win  be  diatribifted. 

Rom  each  dehqr-time  set  die  cone^ionding  reedved  power  can  be  calculated  widi  (7.2)  fcr  any 
frequency.  The  result  is  a  power  value  g  at  frequency  f  with  probability  density  p.^.  To  obtain 

die  power  profadMlity  function  P|(^  at  constam  frequency,  we  should  add  an  paotialnlities  PtOl)  of 

the  tune  sets  i  for  which  the  power  g«  S(f)has  die  same  value: 
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Pg(a)  “  2  Pt(t)  (7.6) 

-  g) 

With  this  method  we  can  calculate  the  distribution  of  the  power  qtectrum.  For  more  detail  on  this 
point  seen  Appendix  C. 

With  the  help  of  a  conqmter  three  typical  values  of  die  power  spectrum  are  calculated; 

1.  the  90%  value;  this  is  the  value  of  S(f)  below  which  90%  of  all  calculated  powers  All, 

2.  the  mean  value;  this  is  the  50%  value, 

3.  the  10%  value. 

hi  the  figures  7.1a..c  the  three  typical  values  of  the  power  spectrum  are  givra  for  the  following 
situation: 

1.  Number  of  paths  is  2 

2.  Ai  =  1;A2  =  0.842 

3.  tj  =  0  )is;  =  0.972  ps 

4.  Delay-time  standard  deviation  o  6.38  ns 
The  carrier  frequency  is  66  MHz. 
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Fig.  7.1c:  ftequeiKy  lesponse  65  -  67.8  MHz 

From  the  figures  7.1a..c  it  can  be  seen  that  at  the  higher  fiteqoencies  tfaeie  is  a  substantial  gap 
between  the  0.10  and  0.90  fraction  curves.  At  lower  frequencies,  fliis  difference  becomes  smaller. 
The  g^  depends  on  die  frequency,  since  at  higher  frequencies  die  delay-time  distiibutioa  causes 
greater  differences  in  the  argument  die  sine  and  cosine  functions  of  (7J2),  whereas  at  lower 
fiequenciet  the  del^-time  uncertainty  causes  a  relatively  small  phase  vasiatkm. 

At  a  frequency  of  70  MHz,  a  full  period  of  2x  is  reached  already  by  a  debqr-time  change  of  At  > 
14.3  ns.  Since  this  time  is  about  twice  the  deli^-time  standard  deviation,  an  possible  power  values 
between  (Aj  -  A2^  >  -16  dB  and  (A|  +  *  5.3  dB  win  occur.  Consequently  there  win  be  a 

great  spread  in  die  resulting  power  distiibutioa  at  fiequcndes.  At  die  relatively  low 
fiequency  of  30  MHz  diis  is  not  die  case,  dierefoie  the  resulting  pow«r  values  win  be  more 
concentrated. 

Hiis  method  shows  us  that  die  measured  ddqr-time  differences  are  not  accurate  enough  for 
calmlating  die  relalive  RF-phase  relttkms  between  the  different  padis  over  a  wdde  frequency 
range. 
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However,  die  measured  RF-phase  values  between  the  different  padis  are  known  at  the 
measurement  frequency.  This  infrMmation  can  be  used  to  restrict  the  delay-time  of  padi  i  to  a 
reladvety  «m«ll  immber  of  discrete  values.  The  RF-|diaae  at  frequency  f  for  path  i  is: 

«(f)  -  -  2afti  (7.7) 

Here  is  the  measured  RF-idiase,  whidi  is  the  RF-phase  at  die  measurement  frequency,  so  tj  is 
restr:''ted  by: 


-  0  (7.8) 

Now  the  probability  density  function  of  each  delay-time  becomes  a  discrete  function.  By  a  non¬ 
linear  transformation  the  time  variable  is  changed  towards  die  nearest  delay-times  x(  diat 
fulfills  (7.8). 

The  time  discretization  is  madiematically  described  by: 


(7.9) 

With  [z]  is  the  integer  so  that  lx-[x]l  is  mmimim  and  ^  is  the  carrier  frequency  of  the 
measuiemenL 

The  probability  of  the  resulting  del^-dmes  is  found  by  additioo  of  aH  probability  densities  of 
die  times  diat  lesuh  in  xn'.  Rgure  72  shows  both  the  analogue  and  discrete  deliy-time  distribution 
in  case  the  radio  channel  consists  of  two  padis. 
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Fig.  7.2:  Delay-tmie  diatributioii  function 

Now  Ae  cotiected  delay-tune  will  result  in  die  same  RF-idiase  ^  as  measured  for  all  times  tj. 
Hence  for  this  fiequency  the  power  density  function  cmsists  of  only  one  power  conqwnent  widi 
value  1.  For  all  other  fiequencies  diis  is  not  the  case.  For  different  delay-times  tj  die  conected 
delay-times  result  in  a  different  power  and  again  a  power  distribution  will  result,  however,  now 
it  will  be  much  less  qiread  especially  around  die  carrier  fiequency. 

The  same  calculations  as  performed  in  figure  7.1  are  now  rqieated  for  die  case  time  discretization 
is  qiplied.  The  results  are  shown  in  figures  73a..c. 
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Fig.  7.3a:  Rrequencyieqxnise  30 -70  MHz.  with  time  diaaetuation 


Pig.  7.3b:  Rnequcocy  n^ooM  30  •  32J  MHz,  wtt  time  diienliaMiaB 
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w  M  •«  *1  n.t  t>* 


Hg.  7.3c:  frequency  leqxmae  65  •  67.8  MHz,  with  thne  diacRlization 

These  figures  show  that  the  results  are  ahnoct  die  same  at  relatively  low  fiequeades.  At  hi|^ 
fiequencies,  however,  the  power  distiftotiaa  is  cooidetefy  dilfaeat  because  of  die  time 
discretization.  At  die  measurement  fieqnency  of  >  66  MHz  dwte  is  no  spread  In  the  distifiioiion 
as  described  in  die  foregoing.  Stanmg  from  die  measnrement  frequency,  die  power  disliibution 
becomes  more  spread. 

7.4  Condnaion 

Using  the  measured  pedi  penuneters,  the  power  speciinm  of  die  frequency  resposiee  can  be 
palrwleled.  Tbe  lesnlla  are  vety  accnrate  when  (be  measured  debqr-lbiies  are  dighlly  cbauged  so 
that  they  are  mneistmt  with  die  measured  relative  RFphaaee.  Thcee  operations  lemlta  in  a 
lunnhw  of  diacrele  delay-time  values  for  every  path.  | 

The  calculated  power  spwctnim  a  quite  aocnraie  over  a  mndi  wider  bendwMfc  than  (he  ! 

bandwidlb  need  fia  die  measmemems,  however  dds  is  dependent  on  fte  dday-dnie  accuracy. 

I, 

I 
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RESULTS  m  A  SIMULATE)  MULTIPATH  ENVRONMENT 


b  ttiis  chapter  dte  remits  of  die  measnremeat  system  operstiiig  in  a  oontndUble  sitnnlsleil 
mnltqMth  eavnaament  are  diaciisaed. 

To  detennine  dte  perfonnance  of  the  measuremeot  qrstem  we  have  simnlated  a  nmlt^padi 
enviroamenL  For  a  maxirmnn  nmidw  of  three  padia  different  snylilndns  and  deb^-times  can  be 
mannally  adjnsted.  Also  it  is  possible  to  add  noise  and  odier  interference  signals.  For  different 
confignrstioiis  and  at  different  frequencies  measurements  are  performed.  t^x»  data  processing 
the  results  ate  compared  widi  die  chosen  parameters. 

8.1  Set-iq>  of  the  simulated  multqiadi  environment 

The  block  diagram  of  die  simnlated  mnltipath  environmerit  is  given  in  fignre  8.1. 


De 

_ 

Pig.  8.1:  Block-dfegiam  of  die  shnulaiednialtfeodi  radio  dmael 


To  simnUte  a  path  of  ttie  mnUqMth  radio  diaimel  we  ue  an  allfflnatnr  and  a  coaxial  cable  of 
certain  length.  The  length  of  die  cable  determines  die  RF^phaae  and  die  delay-time.  The  cable 
lengdi  can  vaiy  iqt  to  about  300  m  (•  U  pa).  The  exact  delay-time  of  a  cable  ia  determined  widun 
an  accnracy  of  -f/-  6  na.  Delay-timea  vaiying  fimn  0  to  1.5  |ia  can  be  lealiaed  in  Ifaia  wiy.  The 
an^ditnde  of  eadi  padi  ia  determined  widiin  -f/- 1  dB  ranging  from  0  to  -12  dB. 

The  dififerent  padi  reqxmaea  are  combined  and  additinnally  noiae  or  an  mterference  aignal  can  be 
added.  The  noiae  power  aa  well  aa  die  power  and  fietjoency  of  die  jammer  can  be  vaiied.  AH 
in^tedances  are  matched  to  avoid  reflectiona  in  one  and  the  aame  cable. 

8.2  MeaaareiiKnta  and  leaulta 

The  meaanrementa  are  divided  into  fonr  categories. 

1.  Accuracy 

2.  Addition  of  an  interfering  signal 

3.  Addition  of  noiae 

4.  Measurement  and  calculation  of  die  fieqnency  reqxmae 

8.2.1  Accuracy 

For  die  first  group  of  measurements  the  system  is  tested  to  determine  the  accuracy  widi  wfaidi 
amplitude  and  del^t-dmea  are  estimated  The  noiae  power  as  well  as  die  jammer  power  are  xero. 
For  equal  configurations  die  measurements  are  executed  at  different  frequenciet  qireod  over  the 
VHF  band.  Fbr  one  fieqnency  die  same  measurement  is  executed  a  number  of  thnea  (13)  to 
examine  the  infhienof  of  time  dqiendent  parameters  of  die  system.  RbsuIib  of  the 
can  be  found  in  the  taUea  A.l  -  A.11  of  ^pendix  A. 

I¥om  these  resulta  it  follows  dmt  die  snqditndes  relative  to  that  of  tte  strongest  path,  esn  be 
determined  addda-f/- 2  dB  for  moneaiea.TI»keM«eipnwil«wWi«m«rim«l«ii|iMHMte»MTnrrwi« 
of  lAi,estimatedlflAi,ootiecti  *  detay-limes  of  the  padis  are  ""ti— <t  widi  an 

accuracy  of  2  sample  instants  to  eidier  side:  Ifitl  S2T,-31.2Sns. 

8.2.2  Additiott  of  an  hMerfateg  signal 

For  the  second  group  of  mranwrinents,  the  inflnence  of  interference  is  examined.  A  jantmer  at  a 
fieqnency  separated  200  kHx  fium  die  canrier  fieqnency  is  used  to  dfasmfc  dethed  signaL  Hk 
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jignal  power  of  Ihe  jaaner  ligiiil  la  of  the  lame  order  of  magnitnde  as  die  desired  aignaL  Results 
of  the  meaaurementi  a«  depicted  ia  table  A.12  of  1^*°**“  ^ 

Rrom  Ifaeae  leanltf  we  can  cooclnde  diat  die  mflamre  of  die  jammer  can  be  neglected  aince  it 
doea  not  canae  atrong  deviationa  from  the  configniation  with  no  jammer  pieaeot 

8.2.3  Addition  of  noiae 

If  we  add  noiae  the  aame  reanlt  ia  found  as  for  addition  of  an  interftring  dgnal.  Low  iqmt-SNR 
(down  to  0  dB)  do  not  infhieoce  die  meaauiement  leanlta  aignificandy.  Reaulta  are  dqncted  in 
table  A.13  of  appendix  A. 

8.2.4  Conversion  of  die  reanlta  to  die  frequency  reqicnae 

The  third  multqiadi  aimulation  inclndea  the  translation  of  the  rnnhipaih  diaracteristica  to  die 
firequen<7  leqxmse,  as  deacrdied  in  duqicer  7.  Fbr  dna  meaanianent  a  configuration  of  three  padis 
is  simulated.  Beside  the  direct  path,  two  extra  padia  were  used: 

Pathl;ti»  Ons.A|*  OdB, 

Path2:T2«  323iis,A2w  OdB, 

Path  3;  T3  «  993  ns.  A3  -  -3.6  <ffi. 

Three  measurementt  were  performed,  eadi  with  a  difierent  padi  configuration:  Padia  1  and  2, 
paths  1  and  3  and  the  configmtion  widi  all  three  padia.  Hxan  die  riimihtiri  algorithm  die 
following  leault  ia  detennmed: 

Pathlixjw  Oiis,Aj«i  OdB,f|B0p, 

Path2:x2B  S16na,A2«  -1  dB,f2«-77®, 

Path 3: X3- 1000 ns, A3-  -6dB,f3--15sO. 


Pkulta  of  die  dday-time  are  doe  to  time  quantization.  A  standard  deviation  of  Oj  -  -  03  -  18 

ns  ia  used.  This  oorremunda  with  d»  ailuation  in  wliicfa  eaiA  of  the  dday-n— -  jg  diatramted 

untform^  over  an  fanerval  of  3  times  die  sample  time  T,  (3T,  -  47  as).  Thon^  foe  meatmementa 

were  executed  on  the  eame  foeqnenctea  m  med  in  — n-i 

an  deecrfoed  only  far  the  caae  ^  -  53.2  MHz.  Fbr  all  ofoer  ftaquenciea.  die  same  madia  me 

in—MB 
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Hg.  8^:  Frequency  response  of  the  radio  channel  made  of  padt  1  and  padi  3,  estimatian 


Flg.8.6:  fteqMacyrefpoBaeofilieraJ»diaHBelinadeofpndil,paat2,aadpaBi3.reiBto 
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Fig.  8.7:  Frequency  response  of  tbe  radio  chaimelnUMile  of  p«&l,  path  2  and  pa&  3,  estimation 
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CONCLUSIONS 

To  determine  the  multipath  characteristics  of  a  radio  channel,  a  measaremeDt  system  has  been 
developed  and  implemented.  This  inclndes  a  transmitter-receiver  configuratioD  using  Direct 
Sequence  Spread  Spectrum  modulation.  The  firequency  band  of  operation  is  dte  VHP-band,  30  - 
70  MHz.  Upon  a  measurement  session,  results  are  stored  on  disc.  Software  data  processing  leads 
to  an  estimate  of  the  nmltipath  characteristics. 

From  the  theoretical  evaluation  and  the  measurement  results  obtained  ftom  a  simulated  mult4>ath 
environment  it  follows  that  the  multipath  characteristics  of  a  radio  channel  are  conqpletely 
described  by  the  number  of  paths  and  their  respective  time-,  phase-  and  delay-time  differences. 

Further  it  turns  out  that  the  technique  presented,  using  a  PN  nK>dulated  carrier  signal  which  is 
coherently  detected  at  the  receiver,  serves  best  the  goal  to  obtain  the  delay-times  as  well  as  die 
conqilex  amplitudes  of  aU  paths  in  a  multipath  environment. 

By  using  PN  sequences  the  demands  on  power  requirements  of  the  transmitter  are  limited. 
Further,  by  processing  a  number  of  consecutive  PN  sequences  it  is  feasible  to  increase  the  sanqile 
rate  virtually  by  a  factor  equal  to  this  number  of  sequences. 

By  applying  a  double  correlatioa  procedure,  accurate  results  are  obtained.  At  the  first  correlation 
stage,  the  received  data  is  correlated  widi  the  PN  sequence  itself.  In  the  second  processing  stage 
the  obtained  correlation  function  is  correlated  with  the  reference  function. 

The  errors  for  each  path  due  to  noise  for  an  iiqnit  signal-to-noise  ratio  of  SNR  s  -14  dB,  are 
theoretically  found  to  be: 
delay-time  :  22  ns 
amplitude  :  0.5  dB 
RF-phase  ;  3  degree; 

Estimation  of  the  multipath  characteristics  of  a  simulated  multqiath  radio  channel  widi  known 
path  delays  and  artqilitude  losses  results  for  iqnit-SNR  »  0  dB  in  the  following  etiurs: 
delay-time  :  ±  7.8125  ns 
amplitude  ;i2dB, 
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The  mult4)ath  characteristics  of  the  radio  channel  can  be  converted  to  die  fireqnency  dtunain.  The 
narrowband  frequency  response  of  the  radio  channel  is  found  widiin  several  dB'g  around  the 
measureinent  frequency.  Theoretical  and  measurenaent  results  of  die  narrowband  frequency 
response  of  a  simulaled  three  path  radio  channel  agree  widiin  ±1  dB  for  a  frequency  band  of  10 
MHz  centred  around  the  rneasurement  frequency. 
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ABBREVIATIONS 

A/D  :  Analogue-Digital 

A^  :  Canier  an^litnde 

1^  :  Contents  of  stage  k  of  shift  register 

Aj  :  An^litude  path  i 

b|^  :  Symbol  k  of  sequence 

Bj^  :  Noise  Bandwidth 

D  :  Number  of  sequences 

Ep  :  Pulse  Energy 

f(  :  Transmitter  carrier  frequency 

h(t)  :  Impulse  response 

I  :  In-phase 

K  :  Sequence  length 

L  :  Shift  register  length 

M  :  Number  of  saoqiles  per  symbol 

N  :  Number  of  paths 

Nq  :  Noise  spectral  density 

P  :  Number  of  sanqtles  under  consideration 

PLL  :  Phase  Locked  Loop 

PN  ;  Pseudo-Noise 

Q  :  Quadrature  phase 

R  :  Correlatiao  function 

R(t)  :  Continuous  conelation  function 

B(m)  ;  Discrete  cooelation  functimi 

r  :  Received  signal 

R^  rCh^ratel/rg 

RF  :  Radio  frequency 

R,  :Sanq>]entel/r, 

S  :  Signalpower 

s  :  Ttansmitted  signal 

SNR  ;  Signal  to  Noise  Ratio 

T  :  Symbol  duration 

t  iTime 
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Tf  :  Chip  duration 

Tg  :  San^)ljng  time 

VHF  :  Veiy  High  Frequency 

Error  value 
Phase  shift  of  path  i 
Variance 

Delay-time  of  path  i 
Radian  fiequency 


♦i 

<?■ 

■h 


(group  leader) 


Ir.  GJM.  Janssen 
(autiior) 
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TABLES  OP  RESULTS  FROM  MEASUREMENTS  IN  HIE  SIMULATED  MULTIPATH 
ENVIRONMENT 

Configuration: 

2padis 

Padi  1;  tj  s  0  8  A}  s  0  dB 

Path  2;  t2  *  S23  ns  A2  *  0  dB 
SNR«40dB 


Table  A.1:  Measurement  results 


Frequency  Delay-times  (ns)  Amplitude  (dB) 


(MHz) 

Conect 

Measured 

Conect 

Measured 

30.425 

523 

531 

0 

■4 

39.975 

523 

516 

0 

-3 

44.475 

523 

500 

0 

4l 

49.825 

523 

531 

0 

-4 

53.200 

523 

516 

0 

-1 

70.925 

523 

516 

0 

-2 
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Cttifigiiration; 

2paths 

Padil:t]>0s  A|-OdB 

Padi  2;  t2  =  993  na  A2  =  0  dB 
SNR-40  dB 


Table  A.2;  Measnrement  results 


Frequency  Dehty-dmes  (ns)  Anq>Iitude  (dB) 


(MHz) 

Correct 

Measured 

Correct 

Measured 

30.425 

993 

1016 

0 

-3 

39.975 

993 

1016 

0 

+0 

44.475 

993 

969 

0 

40 

49.825 

993 

1016 

0 

-1 

53.200 

993 

1016 

0 

+2 

70.925 

993 

938 

0 

42 
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Coofigundon: 

Spadis 

Pitbl:t^BOs  A^oOdB 

Pad)  2:  IS  16  ns  A2~0dB 
SNR-40  dB 


Table  A.3:  Measurement  results 

i^reqiiency 

Deli^-times  (ns) 

Amplitude  (dB) 

(MHz) 

Correct 

Measured 

Correct 

Measured 

30.42S 

1516 

1484 

0 

-3 

39.975 

1516 

1531 

0 

-2 

44.475 

1516 

1516 

0 

+0 

49.825 

1516 

1516 

0 

-6 

53.200 

1516 

1516 

0 

+1 

70.925 

1516 

1516 

0 

-4 

I 


fi 

f 


r 
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Configuratioa: 

2patbs 

Paiiil:ti>Os  Aj-OdB 

Padi  2:  t2  *  523  ns  A2  *  -4  dB 
SNR«40dB 


Table  A.4:  Measmement  results 

Frequency 

Delay-tunes  (ns) 

Amplitude  (dB) 

(MHz) 

Correct 

Measured 

COtrect 

Measured 

30.425 

523 

500 

-4 

-2 

39.975 

523 

547 

-4 

-4 

44.475 

523 

531 

-4 

+1 

49.825 

523 

516 

■A 

-2 

53.200 

523 

516 

■4 

-2 

70.925 

523 

516 

-4 

■6 

I 


1 
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C(nfiguration: 

2padis 

Pathl:tjsOs  Ai=OdB 

Path  2:  t2»  993  ns  A2  =  -4dB 
SNR-40  dB 


Table  A.S:  Measurement  results 


Piequency 

Delay-times  (ns) 

Amplitude  (dB) 

(MHz) 

Correct 

Measured 

Coaect 

Measured 

30.425 

993 

1000 

■4 

-5 

39.975 

993 

1016 

■4 

-3 

44.475 

993 

984 

■4 

-3 

49.825 

993 

1000 

■4 

-1 

53.200 

993 

1000 

-4 

-6 

70.925 

993 

1000 

■4 

-6 

lii 
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Configuration: 

2padis 

l:tj*0  8  Aj=OdB 

2:t2  =  1316tts  A2S-4<1B 

x40dB 


Table  A.6:  Measurement  results 

Frequency 

Delay-times  (ns) 

Amplitude  (dB) 

(MHz) 

Conrect 

Measured 

C^ortect 

Measured 

30.425 

1516 

1531 

-4 

-3 

39.975 

1516 

1516 

-4 

-3 

44.475 

1516 

1500 

-4 

-5 

49.825 

1516 

1516 

-4 

-4 

53.200 

1516 

1531 

-4 

-4 

70.925 

1516 

1531 

-4 

-5 

i! 

I 
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Configuratioo: 

2padi8 

Pathl:t^KOs  Aj«OdB 

Padi  2:  t2  «  523  ns  A2  =  -10  dB 
SNR  =  40dB 


Table  A.7;  Measurement  results 


Frequency 

(MHz) 

Delay-times  (ns) 

Correct  Measured 

Amplitude  (dB) 

Conect  Measured 

30.425 

523 

547 

-10 

-12 

39.975 

523 

516 

-10 

-6 

44.475 

523 

516 

-10 

-16 

49.825 

523 

547 

-10 

-12 

53.200 

523 

500 

-10 

-11 

70.925 

523 

531 

-10 

-10 

TNOraport 


Appendix  A 


Page 

A.8 


Configuration; 

2  paths 

Pathl;tjsOs  Ax  =  OdB 

Path  2;  t2  =  993  ns  A2  =  -10  dB 
SNR  =  40dB 


Table  A.8:  Measurement  results 


Fiequency  Delay-times  (ns)  An^Utude  (dB) 


(MHz) 

Correct 

Measured 

Conect 

Measured 

30.425 

993 

969 

-10 

-10 

39.975 

993 

1000 

-10 

-12 

44.475 

993 

953 

-10 

-8 

49.825 

993 

938 

-10 

-12 

53.200 

993 

953 

-10 

-9 

70.925 

993 

1031 

-10 

-8 
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Configuration: 

2  paths 

Pathl:tj=0$  Aj-OdB 

Path  2:  t2=  IS  16  ns  A2  =  -10dB 
SNRr=40dB 


Table  A.9:  Measurement  results 


Frequency 

(MHz) 

Delay-titnes  (ns) 

Correct  Measured 

Arqtlitude  (dB) 

Correct  Measured 

30.425 

1516 

1500 

-10 

-13 

39.975 

1516 

1516 

-10 

-12 

44.475 

1516 

1484 

-10 

-10 

49.825 

1516 

1531 

-10 

-11 

53.200 

1516 

1500 

-10 

-12 

70.925 

1516 

1547 

-10 

-12 
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Configuration: 

2  paths 

Pathl:t^=Os  Aj  =  0dB 

Path  2:  t2=  1516  ns  A2  =  -10dB 
SNR  =  40dB 


Table  A.  10:  Measurement  results 


Ftequency  Delay-times  (ns)  Anq>Utude  (dB) 

(MHz)  Correct  Measured  Conect  Measured 


53.200 

1516 

1516 

-10 

-13 

53.200 

1516 

1516 

-10 

-13 

53.200 

1516 

1531 

-10 

-13 

53.200 

1516 

1531 

-10 

-11 

53.200 

1516 

1531 

-10 

-11 

53.200 

1516 

1531 

-10 

-11 

53.200 

1516 

1531 

-10 

-11 

53.200 

1516 

1531 

-10 

-11 

53.200 

1516 

1500 

-10 

-11 

53.200 

1516 

1500 

-10 

-11 

53.200 

1516 

1516 

-10 

-11 

53.200 

1516 

1500 

-10 

-11 

53.200 

1516 

1516 

-10 

-11 
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Configuration: 

3  padis 

Path  1:  tj  =  0s 
Path  2:  t2  =  S23  ns 
Path  3:  tj  -  993  ns 
SNR»40dB 


1 


Ai  =  0dB 
A2  =  0  dB 
A3  =  -5.6dB 


Table  A.  1 1 :  Measurement  lesttlts 


Frequency 

Delay-tunes  (ns) 

Anq>litude  (dB) 

(MHz) 

Cotrect 

Measured 

Cotrect 

Measured 

30.425 

523 

516 

0 

1 

-2 

993 

938 

-5.6 

-8 

39.975 

523 

516 

0 

-3  i 

993 

1016 

-5.6 

-6 

44.475 

523 

516 

0 

+3 

993 

1016 

-5.6 

-6 

49.825 

523 

531 

0 

-2 

993 

969 

-5.6 

-8 

53.200 

523 

516 

0 

-1 

993 

969 

-5.6 

-6  ' 

70.925 

523 

547 

0 

-2 

993 

984 

-5.6 

-7  ’ 

I 
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Configuratioo; 

2  paths 

Padil:tj=Os  Aj=OdB 

Path  2:  t2=  1516  ns  A2  =  0dB 


^jammer  =  ^^-^25  MHz 

SNRjaniina'  ”  ^PN^jammer  =  ■•  <1® 
SNR  K  40  dB  (no  jammer) 


Table  A.12:  Measurement  results 


Fluency  SNRj,n,tnw  Delay-times  (ns)  An^litude  (dB) 
(MHz)  (dB)  Correct  Measured  Correct  Measured 


49.825 

15 

1516 

1531 

0 

-2 

49.825 

10 

1516 

1531 

0 

-3 

49.825 

5 

1516 

1500 

0 

+15 

49.825 

0 

1516 

1500 

0 

-50 
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Configuration: 

2paths 

Paliil:tj=Os  Aj  =  0dB 

Path  2:  t2=  IS  16  ns  A2^0dB 
SNR  =  ..dB 


Tabie  A.13:  Measurement  results 


Frequency 

(MHz) 

jammer 

(dB) 

Delay-times  (ns) 

Correct  Measured 

Amplitude  (dB) 

Conect  Measured 

49.825 

5 

1516 

1500 

0 

-5 

49.825 

5 

1516 

1500 

0 

-4 

49.825 

5 

1516 

1500 

0 

-3 

49.825 

5 

1516 

1500 

0 

-2 

49.825 

5 

1516 

1500 

0 

-2 

49.825 

5 

1516 

1531 

0 

-2 

49.825 

5 

1516 

1500 

0 

-4 

49.825 

5 

1516 

1484 

0 

-4 

49.825 

5 

1516 

1531 

0 

-3 

49.825 

5 

1516 

1484 

0 

-1 

49.825 

0 

1516 

1500 

0 

-3 

49.825 

0 

1516 

1500 

0 

-3 

49.825 

0 

1516 

1484 

0 

-6 

49.825 

0 

1516 

1484 

0 

-6 

49.825 

0 

1516 

1500 

0 

-4 

49.825 

0 

1516 

1500 

0 

-3 

49.825 

0 

1516 

1484 

0 

-3 

49.825 

0 

1516 

1500 

0 

-3 

49.825 

0 

1516 

1484 

0 

-4 
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CALCULATION  OF  THE  MULTIPATH  CHARACTERISTICS 


1 

Delay-dme  detennination 

B^ 

2 

Parameter  optiimzati<ni 

B3 

3 

Con^Iexity  and  convergence  of  the  algoridim 

B.9 

'RK>rapMt 
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In  (S.3'j  and  (SJ)  die  ex{BV8aiooa  (4.31)  and  (4.32)  for  ipe]  and  (}[k]  leqiectivdy  can  be 
aubstinited.  By  interdianging  padi  addilion  and  conelatim  we  obiain  the  following  eiqneasioos: 

••  H 

Rlinl  -  ^  I  {AciF«ii(In+k]TM5-^i)  + 

AaiPllQi(tn+k]T*D-^i>  +  "j  ( (n+kJTj^)  )xtkj  - 

N  ~ 

-  S  Z  {AciP»ii(ln+k)T*c-Ti)  + 

‘aif^Qlltn+kJ^AD-^i)  +  "j  ( tn+klTj^,)  )x[k]  - 

-  I  Z  (AciSKiidn+k-lilTju,)  + 

i-l  k— • 

“  *  >*i  tn)  p  (B.l) 

K 

RQ[n]  -  ^Z^-A,iRgQ[n-lil  +  Aj-iRigln-lj^]  +nQ[n].  (B.2) 

In  (B.l)  and  (B.2)  the  following  definitiona  are  used; 


Bjl [n] 

-^_Z  PHii([n+kIT,U,)x[k), 

(B.3) 

RQl[n] 

-^_^PHqi  ( tn+kJT,^^)x[k] , 

(B.4) 

-^_Z_PHiq  <  tn+kj  Tj^,)  X  Ik) , 

(B.5) 

BqqI**] 

-  Z  PHnndn+kJI^lxlk), 

k— -«• 

(B.«) 
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nifn]  -  r  nj{[n+klTj^ij)x[kI, 

(B.7) 

m 

ngln]  X  ngt  [n-fk]Tju))x[k], 

(B.8) 

li  -  int(ti/Tju,) 

(B.9) 

In  (B.9)  int(x)  denotes  the  integer  so  dial  0x-int(x)0  is  miniitnnn.  Widi  dtis  fotnmla  the  analogue 
time  T|  is  converted  into  die  disci^  time  variable  Ij.  Note  that  the  correlation  functions  Rji[n], 
RQ|[n],  Rqiin]  and  RQQ[n]  (which  are  called  die  reference  correlation  functions  in  duqiter  5) 
are  determined  only  by  the  system  and  are  not  dependent  tm  the  radio  channel  c.q.  mnirtpurti 
characteristics  itself.  The  correlation  functions  RqiM  and  RiQ[n]  are  due  to  cross-talk  between 
the  I-  and  Q-channels.  These  functions  are  used  as  reference  functions  to  estimate  die  maltqiadi 
characteristics.  These  reference  functions  are  measured  under  laboratory  conditions  with  die 
transmitter  and  receiver  connected  diiectiy  together. 

I^om  (B.l)  and  (B2)  ve  can  see  diat  the  calculated  correlation  function  is  determined  by 
the  multqMdi  ebaractmixics  and  Ij  for  i  >  1^,..  JI; 

the  noise  terms  nj(n]  and  ngln]; 
the  system  via  RQ[n],  Rgitn],  R^[n]  and  RqqM. 

We  can  now  define  die  criterion  for  die  data  processing  to  Mtimte-  the  mnltipitii  characteristics: 
Determine  the  number  of  paths  (N)  and  fat  each  padi  die  in  die  I-  and  (^-channel  A^, 

A|i  and  delay-time  Ij  so  diat  the  energy  contained  in  all  ratmatM  of  die  noise  sangiles  n][n]  and 
ogCn]  is  minimom.  In  odier  words;  find  N,  A^j.  A^  and  Ij  for  i^l^viN  so  diat 

a2  -  I  AIn]2  -  Aj*  +  Aq2  -  I  (Ailn])*  +  I  (AQ[n])2  - 
n— «  !)•-«  n— • 

*  ■  j|j**ai*QQl“-li]  +  AciRiQ[n-lil  )*,  (B.IO) 

is  miniinunt.  Since  the  oieigy  of  i^[n]  mad  RQ[n]  is  ooocentraled  in  die  pe^  of  die  fimctiaB 
only,  the  WHWiyuioii  of  (B.IO)  can  be  resliicled  to  a  fimited  nanber  of  P  sam|des  starthig  at  p. 
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Tbe  panuneten  p  and  P  must  be  detenmned  fiist  and  differ  for  each  radio  channeL  Any  dioice  in 
which  the  P  selected  samples  cover  most  energy  of  foe  coiielation  fonetkm  satisfies  but  since  foe 
processing  conqilexity  linearly  dq^kls  on  P,  the  parameter  should  be  chosen  as  small  as 
possible. 

An  iterating  algmifom  is  used  to  ddeimitte  foe  mnltipafo  characteristics.  The  block  stmetnte  of 
foe  algorithm  is  shown  by  foe  flowchart  in  figure  B.l.  The  algorifom  starts  wifo  N^l  pafos  and 
optimizes  the  amplitudes  A^.]  and  A,]  and  foe  delay-time  1^.  From  foese  parameters  and  the 
coiielatioa  functions  Rq[o],  Riqln],  RQjIn]  and  RQQ[n],  foe  etror  i?  of  (B.IO)  is  calculated.  If 
this  etror  exceeds  a  certain  threshold  foe  energy  of  die  error  function  is  siqiposed  to  originate 
from  another  path  instead  of  the  noise  terms.  In  diis  case  a  following  path  is  added  and  die 
algoridim  starts  over  again  to  find  the  optimum  solution  for  one  extra  path.  Paths  are  added  until 
the  resulting  error  becomes  less  than  the  error  threshold  value  The  most  important  charts 
of  the  flowchart  of  figure  B.l  will  be  described  now. 


Hg.  B.l:  FlowdiaitoffoealgatUim 
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1  Delay-lime  detBmmitiao 

To  ih^  <nKinnim  atwpiiiiiife  eod  delfly-time  of  an  added  path  the  algotilhm  atatta  by 

»«<itn«»iiig  the  delay-time  l|f.  Thia  time  ia  detennined  from  die  enn  ftmctioo  A^[n]: 

A2[n)  -  (  Rjln]  -  >  + 

H-1  , 

+  {  Rflln]  -  ^X^-AgiRQQtn-lil  +  AoiRiQ[n-li)  )Z.  (8.11) 

The  enor  function  A^(n}  dins  foUows  dhectly  bon  the  lecehred  samples  (R]{n]  and  RQ[n])  and 
the  earlier  estimated  patameten  A^j,  Ajj  and  Ij  for  i  >=  1  A-.-N-l.  Rrom  die  emw  fiinctiiia  A2[n) 
the  maximum  A^[m]  ia  determmed; 

A^Cm]  h  A^In]  fox  all  n,  (B.12) 

with  n  e  (p.p+P-1]  and  m  «  ^jiaP-1].  The  delay-time  ^  is  now  chosen  so  that  the  refoenoe 
conelation  funcdon  R|]%n-l|i|]  and  A^n]  have  their  maxima  on  die  same  instant.  Note  diat  Iji^  is  a 
first  esthnatf  of  the  delay-time  and  may  te  changed  by  the  folhming  st^  of  die  algoaidmi. 

2  Parameser  optiniiiathm 

The  pimeter  eptimiiatiwi  mchidet  bding  the  beat  values  for  A^,  A^^  and  Ij  bar  i  •  1X.|N- 
1  Jf.  Nose  that  the  earlier  found  vafaiea  fin  Agj,  A^j  and  Ij  for  i  •  1.2....N-1  tmqr  be  changed.  The 
opcimiaation  oiaaials  of  an  iaeradnt  algosttm  shown  in  figure  B.2.  The  algorithm  hi  prind^ 
contiaH  of  one  step;  the  drtamnnirinii  of  the  beat  time  shift  ^  <*  (AliAl2f.AlN)  idative  to  die 
current  delqr-lime  set  1  -  (1|  of  the  titne  shift  of  the  delay-time  of  one 

padi  can  take  on  one  of  three  valnea;  -Angp.  0  or  t-A,!^  The  value  of  A^  changea  upon 
program  eaecsrtkn  asri  depands  upon  fte  tnahber  of  paths  N.  bhiailly  A^  -  12,  A^  •  10, 
A,!^  a|...  forNol,  N>2,Na3„  pada.  For  dds  time  riiift  all  poaaiUe  dday-time  seta  ate 
finmed  by  iacreasiag  and  decreaaiag  aaA  delay  dement  wlA  Oe  current  defagr  step  Aglip.  Thit 
mauka  hi  n3  dduHima  sate;  Oi-A^ 

Oi-tAgSip.  l2'^aqp. ..  %A,l,p). ..  Ijadgm^  .  %+A^.  For  each  delay-time  set  tte 

tgirtiianii  tswpHludra  Ag|,Ap2f -.A^  and  A,|,A,2>~’A|N  **  t^lt>s1atad  This  ia  acewapBahad  by 
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cettiiig  die  paitial  derivative  of  die  enor  fimcdcMis  A|2  and  Aq2  of  (B.IO)  be  eqnal  to  zoo.  By 
doing  so,  2N  eqattums  moat  be  lohred  to  obtain  tbe  value  of  the  2N  variablea 
««•  For  dieae  values  of  die  amplitude  the  eator  fnnctian  (B.IO)  is 

given  die  cunent  detay-tane  set  1 = (Ij  J2.-4n)- 


Hg.  B.2:  Flowchatt  on  parameter  opdmizatian 


/ 
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The  mattienutical  derivadoo  of  die  foregoing  is  described  now.  The  partial  derivative  of  CB.IO) 
relative  to  A(.j  is  given  by: 

P+p-1  - 

d/dAc-l  ^  {Afn]>2  -  0 

^  n-p 

P+P-1  -  P+£-l 

X  d/3A(..4{Ailnl)2  +  X  3/dAg.»(AQ[n]  -  0 

n-p  •’  n«p  ■' 

p+p-1  H 

2  {X  +  Rj/j(n~lilRjQ(n-la] )  + 

n-p  i-1  j  »  »  j 

N 

Agi{RQj(n-lj^]Rjjtn-lj]  -  RQQtn-l^lRjQ[n-l  j) )  + 

-  RitnlRiiin-lj]  -  RQ[n]RiQCn-lj])  -  0.  (B.13) 

For  Agj  a  similar  result  is  found: 

P+P-1  - 

3/aA,.,  2  (AfnJ>2  -  0 

n«p 

P+P-1  «  P+P-1 

I  d/dA^j  {Ajtn))2  +  Z  d/dAgj  -  0 

n»p  n**p 

P+P-1  » 

”  dQQlP-lilPQQln-lj] )  + 

-  RiinJRQitn-lj]  +  Rg[nl»QQln-ljl )  -  0.  (B.H) 
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We  now  define  the  iiqjroduct  elements: 


*i,  j 

P4-P-1 

*  2  {Rjl  ^  Rjq Rjq In*l j ]  ) 

(B.15) 

n«p 

p+P-1 

X  {Rqi  jl  ) 

(B.16) 

n-p 

p+P-1 

■  2  {Rlltn-lilRQlln*ljl  -  RiQ[n-lilRQQtn-ljl } 

(B.17) 

n-p 

°i,J 

p+P-1 

■  2  <RQi[n“lilSQl[n-lj]  +  R^[n-lilRQQln-lj] } 

(B.18) 

p+P-1 

■  2  {Ritn“li3RiiI*i*lj3  +  RQ[n“li3RiQ[*'“lj3  > 

(B.19) 

n*p 

p+P-1 

-  2  {RiCn-lilP^lln-lj]  -  RQ[n-li]BQQ[n-lj]> 

n*p 

(B.20) 

Fw  (B.13)  and  (B.14)  we  can  write: 

£  lAci  *ij  ^ai  ^  (B.21) 

[A^i  Cij  +  A,i  Dij]  -  fj  for  j  -  1,2,..,H.  (B.22) 

The  equations  (BJl)  and  (B^)  can  be  written  in  malrix  fiosm  fee  an  j  > 


All  *21 

-  *»! 

®11  ®21 

•  •  •  •  *111 

'*el  ' 

*1' 

*12  *22 

....  Aj^ 

®12  ®22 

....  ®IK 

*c2 

*2 

*111  *211 

•  >  ■  ■  ^WH 

®m  ®2H 

....  ®WI 

*^ 

®« 

<^11  ^21 

....  C|,1 

®11  ®21 

....  Dui 

*al 

®1 

Ci2  C22 

....  °1I2 

®12  ®22 

....  D||2 

*a2 

®2 

>  *^1H  ®2tl 

••••  ®WI 

®lll  ®2JI 

-  ®lw; 

I.  *8*  - 
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The  matrix  eqnatk»  (B.23)  is  solved  1^  LU  decMiqiosition.  die  results  of  die  snylitmiM 
Af  1  Ac2— •^cN  ^  function  of  (B.IO)  is  cslcnUted.  For  one  of  die 

delay-time  sets  the  enor  function  is  minimum.  This  result  in  the  best  delay-time  shift  and  die  new 
best  delay-time  set 

In  case  a  different  best  time  set  results,  ^  »  (Al]Al2...Aliij)  *  (0,0,..,0),  the  algorithm  continues 
by  kxddng  for  die  new  best  delay-time  shift  In  case  the  best  deli^-time  shift  set  is  zero  fin  all 
shifts  (A1j,AI2,..A1n)  =  (0,0,..,0)  die  algaridun  continues  dqiending  tm  die  riiifl  stq>  size  If 
the  shift  step  size  is  minimum,  =  I,  the  algoridim  ends  and  the  cutrent  del^-time  set  1  is  the 
optimum  delay-time  set  for  N  padis.  In  case  ^  1 ,  die  shift  step  size  is  decreased  to  be  able  to 
examine  the  inteimediate  delay-time  sets. 

The  last  stq>  to  obtain  all  mult^iath  cbaiactetiatics  is  the  deteimination  of  die  phase  difieroices  of 
die  found  paths.  These  follow  directly  from  (4. 19)  and  (4.20): 


-  aret«n(A,i/Aoi) . 


(B.24) 


3  Conqilexity  and  convergence  of  die  algorithm 

-  Coiiqfilexity 

Because  die  number  of  paths  is  of  order  of  5  for  most  cases,  the  conqilexity  of  die  anqditude 
calculation  is  not  determined  by  solving  die  matrix  (B.23)  but  by  die  in|aodiict  calculation  (B.15) 

-  (B20).  For  each  mfsoduct,  2P  multqilicatioiis  and  P  summations  have  to  be  r«inii«iw<i  ppr  tiM« 
reason  die  consiaut  P,  die  number  of  sanqile  instants  under  copsidention,  should  be  kqit  as  ««wii 
aspoasiUe. 

-Convergence 

The  algorithm  converges  if  the  enor  A^  decreases  s{kni  progiam  execntion.  If  die  enor  is 
smaller  than  die  dneshoid  vahie  die  algorithm  ends.  During  program  execution  die  emr 

A^  dmges  at  two  locationa  in  die  aigocidun.  Pint,  if  a  new  pedi  is  added  and  second  if  die  ddi^- 

time  set  i  «  (liJ2*"<h4)  *  flanged.  Since  dnmges  in  die  ddqr-thne  set  are  deiennined  by  dr 
enor  A^  itadf,  die  error  win  not  incrfase  iqion  optimizatian  of  die  ddqr-time  set 
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The  deb^-dme  of  a  new  extra  path  ia  detennined  by  the  maximum  of  die  enor  fonction  A^[n}  aa 
deaciibed  in  aecdon  B.l.  The  new  path  canaea  a  deoeaae  of  die  aror  if  die  eneigy  in  the  etrar 
function  A^[n]  ordinates  from  a  path  of  the  radio  channeL  In  caae  die  energy  of  die  enor  fimcdon 
is  caused  by  noise,  die  change  of  the  error  A^  is  undetennined.  In  pnaidce  it  turns  out  diat  the 
error  A^  decreases  alwi^s.  Hence  we  may  state  diat  the  algorithm  converges. 
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CALCULATION  OF  THE  CUMULATIVE  DISTRIBUTION  FUNCTION  OF  THE 
FREQUENCY  RESPONSE 

From  each  delay-dme  set  we  can  calculate  die  coneapondittg  received  power  for  any  frequency 
with  (7.2).  The  result  is  a  power  value  g  at  frequency  f  with  probability  density  Pt(^.  To  obtain 
die  power  probability  density  function  pg(g>  at  cmistant  frequency,  we  dionU  add  all  probability 
density  values  Pf(t)  of  the  time  sets  t  for  which  the  power  g  k  S(f)  is  equal: 

Pg(2)  “  I  <C*1) 

(S(£)l^-  g) 

Since  we  cannot  determine  the  inverse  function  of  (7.2)  to  calculate  die  time  set  I  for  a  given 
power  g,  die  power  probability  function  Pg(g)  cannot  be  siiinlified  fhrdier  then  as  given 
according  to  fannnla  (C.l).  However,  we  can  follow  a  statistical  qiproacdi  for  which  we  do  not 
need  diis  inverse  function.  Fdr  instance,  at  frequency  f,  die  mean  power  E(glf]  can  be  calculated 
by: 


E[glf]  -  r*/  S(f)lx,^  Pt<t)  'dtj,.  (C.2) 

— *0  — «o  —  — 

With  the  aid  of  contyuter  calculations,  we  can  ipproximate  (C.l).  To  do  so,  we  quantify  die 
variables  g  and  t.  The  time  intervals  are  chosen  mmH  wmngii  so  ii««t  the  energy  upectnim  S(f)  for 
all  delay-times  widiin  die  intervals  differs  only  sli^tly.  The  time  scales  are  resided  fiodier  only 
for  die  area  widi  substantial  ptobdbQity  density  of  die  delay-time.  If  we  restrict  die  time  seek  of 
each  dehty-time  to  foe  interval  [^-2o,  tf»-2o],  a  fraction  of  0.93  is  covered.  This  scale  is  divided 
imo  M  intervals  of  equal  fcngfo.  The  probability  pj  of  the  delay-time  to  fidl  in  die  interval  j, 
tj  e  [T|-2o44(j-l  io/M,  T|-2(H4(j)niM]  frir  j>l,2,..M  i*  calcnlattd  with  foe  Q-ftmcthm: 
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ti-20+4(  j)0/M 

Pj  “  J  Pt<t)  dt  “ 

Ti-20+4{j-l)0/M 

-  Q(-2+<J-1)/M)  -  Q(-2+j/M),  (C.3) 


with 


Q(x)  -  /  1/V(2k)  0xp(-y^/2)  dy.  (C.4) 

X 

The  Q-fiinction  is  well  tabnlated  or  can  be  calculated  by  function  qtproximaticm  (Abnmowitz,  M. 
and  Stegun,  lA.,  1968).  For  N  padu  this  mediod  leads  to  delay-time  sets.  For  each  delay- 
time  set  the  cooesponding  power  value  (or  power  g)  is  calculated  fiom  (7.2).  The  lesultiiig 
powers  are  tabulated  and  an  qtptoximatioo  of  (C.l)  is  obtained. 
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